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Abstract 


Seasonal fluctuations in the rate of rotation of the Earth are detected when checking 
the best clocks with astronomical time. These variations may be shown to result from 
transfer of momentum between the atmosphere and the globe. As a consequence, the study 
of these fluctuations may bring interesting information in the field of meteorology. 


1. Les progrès réalisés dans la construction 
et la comparaison des garde-temps ont permis 
de déceler des fluctuations relatives öw/o de la 
rotation w de la Terre, dont l’ordre de grandeur 
est de + 1078. Le caractère irrégulier et la 
périodicité annuelle de ces fluctuations (FINCH 
1950, SCHEIBE ET ADELSBERGER 1950, STOYKO 
1936, 1937, 1949, UHINK 1950) suggérent une 
cause d’origine météorologique (VAN DEN 
DUNGEN et al. 1949). 

Les valeurs normales de ôw/« établies par N. 
StoyKo (1951 a) pour le début et le milieu 
de chaque mois sont représentées à la fig. 1 
(courbe I). Le calcul de ces valeurs normales 
est basé sur la comparaison des déterminations 
d'heure effectuées en six observatoires et 

ortant sur une période de 15 années. Notons 
qu’ainsi N. STOYKO a éliminé, en principe, la 
variation erratique et la variation séculaire de 
la rotation du Globe. En effet, pour chacune 
des 15 années, il a calculé les fluctuations moyen- 
nes dw au début et au milieu des douze mois 
de l’année par rapport à la vitesse annuelle 
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moyenne @,, correspondante (dw =a —«,,). Il 
en résulte que la courbe I de la fig. 1 met 
uniquement en évidence les fluctuations de 
période annuelle; elle montre que la rotation 
de la Terre est accélérée de mars à août et 
retardée d’aoüt à mars. 


2. En l'absence de forces extérieures agissant 
sur le système «Terre — Atmosphère», la som- 
me Mc + Ma des moments cinétiques du 
Globe et de l’Atmosphére, autour de l’axe de 
la Terre, demeure constante et, par suite, à 


: fel OMe 
chaque instant f, toute variation an du 
moment cinétique Mc du Globe est compensée 


IM : 
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traire, du moment cinétique M4 de l’atmos- 
phère; d’où 


par une variation 
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où Ic désigne le moment d’inertie du Globe 
terrestre. Dés lors, une fluctuation de la rota- 
tion de la Terre ne peut provenir que d’un 
échange de quantité de mouvement d'Ouest, 


> F. H. VAN DEN DUNGEN, J. F. COX ET J. VAN MIEGHEM 


-20— —=s 


Ô 
Fig. 1. Valeurs moyennes de ee : 
o 


Courbe I. Valeurs normales au début et au milieu du 
mois de N. STOYKO (1951). 


Courbe II. Valeurs calculées 4 partir des échanges de 
quantité de mouvement d’Ouest entre le Globe et 
l'atmosphère terrestre, dans l’hypothèse d’une compensa- 
tion parfaite de ces échanges à l’échelle du Globe, au 
bout d’une année. (F. H. VAN DEN DUNGEN, J. F. Cox et 
J. VAN MIEGHEM, 1950). 
x W. H. Munk et R. L. MILLER, 1950. 


Courbe III. Valeurs mensuelles normales de la fluctua- 
tion du moment cinétique absolu de l'atmosphère. 


entre le Globe terrestre et l'atmosphère, à 
travers la surface de la Terre. 

Mais un échange de quantité de mouvement 
d’W entre l’atmosphère et le Globe ne peut 
être dû qu’au frottement exercé par la surface 
du Globe sur les courants atmosphériques 
zonaux. Dans cet effet de frottement, il faut 
inclure l'effet provenant des différences de 
pression de part et d’autre des chaînes de mon- 
tagnes orientées du Sud vers le Nord. 

Il est aisé d'établir le bilan (van MiEGHEM 
1950) 


Ma 
ot 


nu Lf (pe — pw) r cos Ô 0 (2) 


- ff Tx COS MP Ep» a 


dans lequel t, désigne la composante zonale 
de la tension due au frottement de la surface 
du Globe 2’ sur la couche d'air adjacente; px 


et py, les pressions prises sur un même paral- 
lèle à l'Est et à l'Ouest des montagnes; a, la 
surface limitée par le profil obtenu en projetant, 
à l’aide des parallèles, le volume occupé par les 
montagnes, sur un plan méridien ; y, la latitude 
et r, le rayon vecteur au centre de la Terre. 

3. Grâce aux tables numériques établies par 
V. P. STARR et par ses collaborateurs (1949), à 
l’aide de données météorologiques de surface 
fort nombreuses, il est possible de calculer les 
deux termes au second membre du bilan (2) 


3 ti) 
et finalement d’en déduire — (VAN DEN 


at 
DuNGEN et al. 1950 a). 


Les tables V et VI du rapport précité de 
V.P. Starr fournissent les valeurs mensuelles 
normales de la quantité de mouvementéchangée 
entre le Globe et l’atmosphère au N du paral- 
lèle de 10° N. Rappelons que le flux de la 
quantité de mouvement dW, à travers la 
surface du Globe, est ascendant là où le vent 
dominant est d’E (latitudes tropicales) et, 
descendant là ot le vent dominant est d W 
(latitudes tempérées). Il en résulte, dans chaque 
hémisphère, un flux de quantité de mouvement 
vers les pôles (STARR 1948). Mais, comme les 
courants d’W sont beaucoup plus forts dans 
l'hémisphère Sud que dans l'hémisphère Nord 
et que, d’autre part, l’asymétrie entre les alizés 
des deux hémisphéres est peu accusée, le flux 
ascendant à la surface du Globe («Terre > 
Atmosphère») aux bafses latitudes de l’hémi- 
sphère Nord, n’est pas compensé, dans l’hémi- 
sphère Nord, par le flux descendant («Atmos- 
phère > Terre») des latitudes tempérées. Une 
partie du flux ascendant à la surface du Globe, 
dans les régions tropicales de l'hémisphère 
Nord, traverse le plan de l'équateur (WIDGER 
1949) et est compensé dans les régions tem- 
pérées de l'hémisphère Sud, où une circulation 
d’W plus intense entraîne, à la surface du 
Globe, un flux descendant de quantité de 
mouvement plus grand que dans l'hémisphère 
Nord. 

De ces considérations, on déduit immédiate- 
ment pour les flux verticaux de quantité de 
mouvement d’W à la surface du Globe et pour 
les flux horizontaux méridiens de la même 
quantité, dans l'atmosphère libre, la réparti- 
tion schématique reproduite à la fig. 2a. En 
partant de ce schéma, on est tenté de dessiner 
les lignes de flux du moment cinétique absolu 
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de l'air, dans un plan méridien, conformément 
au schéma de la fig. 2 b. Cependant il convient 
de noter que ce schéma ne peut être exact que 
sil n'existe pas de circulation moyenne dans 
les plans méridiens ou, ce qui revient au même, 
si, dans l'atmosphère libre, seuls des échanges 
horizontaux de quantités de mouvement sont 
possibles (hypothèse géostrophique). Or, si 
l’on se reporte au profil connu de la circula- 
tion zonale (UNIVERSITY OF CHICAGO, 1947), 
on constate qu'au S du Jet Stream» l'intensité 
du vent zonal suit à peu près la loi de la con- 
servation du moment cinétique absolu, ce qui 
révèle l'existence, aux latitudes subtropicales 
et tropicales, d’une circulation moyenne dans 
le plan méridien. Au N du Jet Stream», par 
contre, l'intensité du vent d’W suit une loi 
toute différente; à savoir, en première approxi- 
mation, la loi de la conservation de la com- 
posante verticale de la rotationnelle absolue 
(RossBy, 1947), ce qui indique la prédominance, 
aux latitudes tempérées, d'échanges horizon- 
taux suivant les méridiens. Bref, aux latitudes 
tempérées, les échanges de quantité de mouve- 
ment entre zones de latitudes différentes sont 
essentiellement des échanges horizontaux, tan- 
dis qu'aux latitudes plus basses, il y a, à la 
fois, échange horizontal et circulation méri- 
dienne moyenne (PALMÉN, 1951). On voit 
ainsi que les lignes de flux du moment cinétique 
absolu de l’atmosphère telles qu’elles ont été 
mises en évidence par V. P. STARR (1948) et 
W. K. WIDGER (1949) ne peuvent exister 
qu'aux latitudes tempérées (fig. 2 b). 

4. En admettant que la compensation entre 
les Aux ascendants et descendants de quantité 
de mouvement d’W, à travers la surface du 
Globe, soit parfaite à l'échelle de la Terre, au 
bout d’une année, nous avons obtenu (VAN DEN 
DuncEN et al., 1950 a), pour la fluctuation de 
la rotation du Globe, des valeurs mensuelles 
moyennes (voir courbe II de la fig. 1) en bon 
accord avec les valeurs déduites des détermina- 
tions de l'heure (STOYKO 1936, 1937, 1949, 
1951 a; SCHEIBE ET ADELSBERGER, FINCH, UHINK 
1950). Remarquons que les amplitudes des 
courbes I et II de la fig. 1 sont du méme ordre 
de grandeur et que, de plus, elles oscillent a 
peu près en phase. Enfin, nous avons complété 
la fig. 1 par deux points (x) qui mettent en 
évidence la variation saisonnière de la rotation 
de la Terre, telle qu’elle a été déduite, par 
W. H. Munk et R. L. Mirter (1950), de la 


Fig. 2. Flux de la quantité de mouvement zonal d’W à 
la surface du Globe et dans l’atmosphère libre. 


différence des moments cinétiques de l’atmos- 
phère en janvier et en juillet. 

En résumé, les fluctuations saisonnières de 
la rotation de la Terre apparaissent comme 
directement liées aux fluctuations des échanges 
de quantité de mouvement d’W entre le Globe 
et l'atmosphère, c’est-à-dire, en dernière ana- 
lyse, aux fluctuations de la circulation atmos- 
phérique zonale. 

5. Inversement, en partant des valeurs nor- 
males des fluctuations de la vitesse de rotation 
de la Terre, établies par N. Sroyxo pour le 
début de chaque mois (StoyKo, 1951 a), on 
peut déduire, grâce à (1), les valeurs men- 
suelles normales des fluctuations du moment 
cinétique absolu de l'atmosphère terrestre 
(voir fig. 1, courbe III). On constate que le 
moment cinétique absolu de l’atmosphere 
augmente d'août à mars (hiver boréal, été 
austral) et diminue de mars à août (été boréal, 
hiver austral), ce qui correspond à une forte 
asymétrie entre les circulations zonales des 
atmosphères des deux hémisphères. 

Ensuite, en combinant les valeurs mensuelles 
normales de la fluctuation du moment cinétique 
absolu de l’atmosphère avec les valeurs men- 
suelles normales des quantités de mouvement 
d’W échangées entre la Terre et l'atmosphère 
au Nord du parallèle de 10° N (fig. 3, courbe 
I), déduites des tables de V. P. Sarr (1949), on 
trouve, grace à (1) et (2), les valeurs mensuelles 
normales des quantités de mouvement d’W 
échangées entre la Terre et l'atmosphère au 


Sud du parallèle de 10° N (fig. 3, courbe Il). 
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Fig. 3. Valeurs mensuelles normales du flux de la quantité 
de mouvement d'Ouest à travers la surface du Globe, au 
Nord du parallèle de 10° N (D et au Sud 
de ce parallèle (I). 


Flux (Terre — Atmosphère»: + 
Flux (Atmosphère — Terre»: — 


On voit ainsi qu'au cours de l’année, le mois 
de juin excepté, le Globe terrestre cède de la 
quantité de mouvement à l'atmosphère au Nord 
du parallèle de 10° N (courbe I, fig. 3) et que, 
par contre, au Sud de ce parallèle, c’est l’atmos- 
phère qui cède constamment de la quantité de 
mouvement à la Terre (courbe II de la fig. 3). 
Le flux «Atmosphère — Terre au Sud du 
parallèle de 10° N est minimum en août et 
septembre lorsque, en moyenne, la vitesse de 
rotation de la Terre est maximum; ce flux est 
maximum en février, lorsque la vitesse de 
rotation de la Terre est minimum (VAN DEN 
DUNGEN et al. 1951). Il en est à peu près de 
même pour le flux «Terre > Atmosphère au 
Nord du parallèle de 10° N (fig. 3, courbe I). 

6. L’asymétrie entre les circulations des deux 
hémisphères, mise en évidence par les courbes 
en escaliers de la fig. 3, permet de conclure 
qu il doit exister, dans l'atmosphère, un flux de 
moment cinétique absolu à travers le plan de 
l'équateur, dirigé de l’hémisphère N vers 
l'hémisphère S et que ce flux doit être maxi- 


mum de novembre à février (été austral, hiver 
boréal) et minimum de juin à septembre 
(été boréal, hiver austral). De juin à septembre, 
il y a à peu près compensation, dans chaque 
hémisphère, entre les flux ascendant et descen- 
dant de quantité de mouvement à travers la 
surface du Globe (van DEN DUNGEN et al., 
1951). Ce résultat est confirmé par les valeurs 
mensuelles normales du flux de moment 
cinétique absolu de l’air à travers la surface 
verticale zonale le long du parallèle de 10° N. 
Ces valeurs normales ont été calculées par 
V. P. Starr (dernière colonne des tables I et 
IV du Rapport précité, 1949). Il est regrettable 
qu’en raison de données aérologiques insuf- 
fisantes, nous ne disposions pas encore des 
valeurs correspondantes pour le flux à travers 
le plan de l'équateur. 


7. Il est aisé de calculer les valeurs mensuelles 
normales des fluctuations du moment cinétique 
absolu de l'atmosphère au Nord (courbe I de 
la fig. 4) et au Sud (courbe II de la fig. 4) du 
parallèle de 10° N, à l’aide des courbes I et 
II de la fig. 3 et des valeurs mensuelles normales 
du flux de moment cinétique absolu à travers 
la latitude de 10° N, établies par V. P. STARR 
et ses collaborateurs (1949). Ces fluctuations 
sont plus accusées dans l'hémisphère Sud que 
dans l'hémisphère Nord; elles sont pratique- 
ment constantes de mars à juin et de septembre 
à décembre dans l'hémisphère Nord et d’octo- 
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Fig. 4. Valeurs mensuelles normales de la fluctuation du 
moment cinétique de l’atmosphère au Nord du parallèle 
de 10° N (I) et au Sud de ce parallèle (II). 
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Eq 
Fevrier 


Eq 
Mai 


Novembre 


Fluctuation du moment cinétique absolu de l'atmosphère; 


—— Flux du moment cinétique absolu à travers la surface du Globe et à travers la surface verticale zo- 


nale le long du parallèle de 10° N; 


Flux résultant à travers la surface du Globe; aux époques de l’année où la vitesse de rotation de la 


Terre est maximum (août—septembre) et minimum (février—mars) et aux époques où cette 
vitesse croît (mai—juin) et decroit (novembre—décembre) le plus rapidement; à la ligne unité: 
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bre à février (été austral) dans l'hémisphère Sud 
(VAN DEN DUNGEN et al. 1951). 

8. Afin de mettre en évidence le lien entre 
les fluctuations du bilan du moment cinétique 
absolu de l’atmosph£re et les fluctuations de la 
rotation de la Terre, nous avons représenté ces 
fluctuations graphiquement 4 la fig. 5, aux 
époques de l’année où la vitesse de rotation de 
la Terre est voisine de ces valeurs extrêmes et 
aux époques où cette vitesse varie le plus: 
vitesse minimum en février, accroissement de 
vitesse le plus fort en mai, vitesse maximum 
en août et diminution de vitesse la plus rapide 
en novembre. 

On constate ainsi 1°) que l'échange de quan- 
tité de mouvement entre les atmosphères des 
deux hémisphères est le plus faible lorsque la 
vitesse de rotation de la Terre est maximum 
(août) et le plus fort lorsque cette vitesse est 
minimum (février), 2°) que l'échange de 
quantité de mouvement entre l'atmosphère et 
la Terre, à l'échelle du Globe, est minime (à 
peu près compensation, à la surface du Globe, 
entre le flux ascendant de l'hémisphère N et le 


flux descendant de l’hémisphère S) lorsque la 
vitesse de rotation de la Terre est voisine de 
ces valeurs extrêmes (en février et en août) et 
que, par contre, cet échange est considérable 
(non compensation entre les mêmes flux) 
lorsque la vitesse de rotation de la Terre varie 
le plus rapidement (en mai et en novembre). 

E. PALMÉN a fait observer que si l’on admet 
le point de vue de la circulation méridienne 
(PALMEN 1951), l’on doit s'attendre à un 
échange maximum du moment cinétique vers 
le Sud, à travers le parallèle de 10° N lorsque 
l'équateur thermique occupe sa position la 
plus septentrionale, c’est-à-dire vers le mois 
d'août ct, à un échange minimum, vers le mois 
de février, lorsque l'équateur thermique occupe 
sa position la plus méridionale. Le fait que l’on 
observe le contraire (voir fig. 5) prouve qu’à 
côté de la circulation méridienne moyenne, 
dont on ne peut nier l'existence (voir fin du 
§ 3), il existe, aux basses latitudes, un échange 
horizontal de quantité de mouvement. 

9. Enfin, à la fig. 6, nous avons reproduit 
les valeurs normales de dw/m au début et au 
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Fig. 6. Variation saisonnière de la rotation de la Terre 
pendant les années 1945—1946, d’après N. STOYKO (1951). 
Courbe I: valeurs normales de Üw/«. 


milieu de chaque mois (SToYKO 1951 a, courbe 
I de la fig. 1) ainsi que les valeurs de dw/w au 
milieu de chacun des mois des années 1945 et 
1946 (StoyKo, 1951 b). La fig. 6 représente les 
anomalies dans la variation annuelle de la rota- 
tion de la Terre au cours de ces deux années, 
notamment: un retard supplémentaire dans la 
rotation de la Terre par rapport 4 la normale 
(courbe I), de septembre à février, et une 


avance supplémentaire, de mars à août. En se 
reportant à la fig. 5, on voit que ces anomalies 
dans la fluctuation annuelle de la rotation de la 
Terre, en 1945 et 1946, peuvent étre expliquées 
par un renforcement de la circulation des 
alizés par rapport à la circulation zonale d’W 
au cours de l’hiver boréal et de l’été austral et 
par un renforcement de la circulation zonale 
d’W par rapport à la circulation des alizés au 
cours de l'été boréal et de l’hiver austral. 

Nous sommes donc en droit de con- 
clure que les fluctuations saisonnières de la 
rotation de la Terre peuvent être considérées 
comme un indice intégré de la circulation 
atmosphérique zonale. Dans l’état actuel des 
moyens techniques, les horloges à quartz ne 
permettent la détermination d’une fluctuation 
dw de la rotation w de la Terre que pour une 
période d’au moins dix jours. Cette fluctuation 
représente alors un indice de circulation 
décadaire global. 

L'étude des fluctuations de la rotation 
moyenne annuelle w,, de la Terre (fluctuations 
séculaires de w) et la recherche d’éventuelles 
causes météorologiques, n’est pas sans impli- 
quer, en plus de l’amélioration des garde- 
temps, des progrès dans la définition de l'échelle 
des temps (cf. VAN DEN D'UNGEN et al., 1950 b). 
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Theory of Development and of Thickness Patterns 


By G. C. McVITTIE (Queen Mary College, London) 


(Manuscript received 6 sept. 1951) 


Abstract 


The equations of motion, of continuity, of vorticity and of heat-transfer are transformed from 
coordinates (x, y, 2) to coordinates (x, y, p) without the use of the hydrostatic equation in the 
vertical. An approximation based on the empirical fact that the isobaric surfaces are slightly 
inclined to the horizontal together with the use of dimensionless variables, gives the hydrostatic 
equation and simplifies the fundamental equations by the rejection of many of their terms. 
Thereafter definitions of certain types of geostrophic and nongeostrophic motions lead, respec- 
tively, to (a) Rossby’s potential-vorticity equation; (b) the development and thickness-patterns 
theory of R. C. Sutcliffe, which is discussed in detail; and (c) J. G. Charney’s (1949) treatment of 
the equivalent barotropic atmosphere. These models of atmospheric motions are compared and 


contrasted. 


I. Introduction 


In recent years, many investigations in 
dynamical meteorology have dealt with the 
vorticity equation, which is a formula giving 
the rate of change “following the motion” 
of the vertical component of vorticity of 
the air. Attention has also been directed to the 
part played by the rate of change “following 
the motion” of the pressure, and to the use 
of upper air charts of “thickness patterns”, 
which exhibit the difference in height between 
two isobaric surfaces above some part of the 
Earth’s surface. Amongst investigations of 
this kind is the theory of development and 
thickness patterns due to R. C. SUTCLIFFE 
(1947, 1950), A. G. Forsdyke and their co- 
workers, papers to which we shall refer 
hereafter as SI and SII. The first paper 
deals with the question of development, the 
main formulae being SI (ro), (11) and 
(14); whilst the second is concerned with 
the use of a formula, SII (3), for the local 
rate of change of the thickness between 
two isobaric surfaces. In view of its im- 


portance in current British meteorological 
research, it has been thought worth while to 
analyze again from first principles the mathe- 
matical basis of this theory. The investigation 
given in the following pages seeks to discover 
how far the formulae referred to above de- 
pend on empirical approximations, on one 
hand, and on a priori definitions of certain 
types of motion, presumed to be possible in 
the atmosphere, on the other. The plan of the 
work is as follows: — In par. 3, the fun- 
damental equations of motion, continuity, 
vorticity and heat-transfer are transformed 
from local rectangular coordinates (x, y, z) to 
coordinates (x, y, p), where p is the pressure, 
without using the assumption that the at- 
mosphere is (vertically) in hydrostatic equilib- 
rium. This part of the analysis is a straight- 
forward application of the general theory of 
transformations given by the author (Mc Vittie 
1948, 1949), papers hereafter referred to as MI 
and MII, and therefore most of the mathe- 
matical details are omitted. In par. 4, dimen- 
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sionless variables are introduced and it is 
shown how the hydrostatic approximation is 
derived from a numerical approximation de- 
pending on the fact that the isobaric surfaces 
are only slightly inclined to the horizontal. 
By the same approximation, many terms in the 
fundamental equations can be neglected and 
these equations cleared of all the parameters 
introduced by the conversion to dimension- 
less variables. The vorticity equation is dealt 
with in par. $ and it is shown that, by means 
of a priori definitions, the Rossby potential- 
vorticity equation and the formulae SI (ro), 
(rt) and (14) can be obtained. The latter are 
more restrictive than the former. In par. 6, 
the thickness formula SII (3) is discussed and a 
criterion for the order of magnitude of its 
various terms suggested. Ks pointed out 
that SII (3) cannot be expected, by itself, to 
throw much light on the dynamics of atmos- 
pheric motions. It plays a part analogous to a 
boundary condition and so must be combined 
with the formulae of par. 5 before dynamical 
conclusions can be reached. Finally in par. 
7 we apply the theory to the equivalent- 
barotropic atmosphere and contrast the results 
with those of Sutcliffe’s method. 

The Coriolis parameters, 

= 20cos gy, 1 = 2-0 sin 9, 

are treated throughout as constants, which, at 
first sight, appears to make the treatment less 
general than SUTCLIFFE’S (1947). But, in fact, 
Sutcliffe makes little use of the variation of | 
with latitude. If it is desired to take account 
of the Earth’s curvature, it is better, in the 
present author’s opinion, to abandon local 
rectangular coordinates and use a system 
appropriate to a curved Earth (McVitTIE 1951). 
The theory of par. 3 onwards could easily be 
re-written in terms of such a coordinate 
system. 

The author wishes to emphasize that he is 
concerned with mathematical theory and not 
with practical forecasting. The numerical 
values of meteorological elements, particularly 
in the upper air, are usually known only 
to a limited degree of accuracy. Now if it 
happens that a relation between them 1s hit 
upon which is useful in practical forecasting, 
this fact by itself does not establish the relation 
theoretically : its usefulness may be contingent 


on the inaccuracy of the data. The theoretician 
is then confronted by the problem: how, if at 
all, can the empirically useful relation be 
derived from the fundamental equations of 
hydrodynamics and thermodynamics? To 
argue that a solution of this problem is super- 
fluous is, in the present author’s opinion, to 
relinquish the hope of understanding meteoro- 
logical processes. 


2. Pressure as vertical coordinate. Notation 


At a point O on the Earth’s surface, whose N. 
latitude is @, the conventional local rectangular 
coordinates are set up. The plane Oxy is 
the tangent plane to the Earth’s surface at O, 
Ox is drawn to the East and Oy to the North, 
whilst Oz coincides with the vertical at O. 
Consider an unsteady motion of the atmosphere 
so that an isobaric surface has the equation 


(2.1) 


where C is a constant. If we solve this equation 
for z, we express this coordinate as a function 
of x, y, t and C. Hence for all values of C we 
can write 


HEC 


22 (api) (2.2) 


in which each value of p is associated with one 
isobaric surface. We regard (2.2) as an equation 
which permits us to replace (x, y, z) by (x, 
y, p) as the coordinates describing the position 
at time f of an elementary mass of air. Moreover 
any function, f, of (x, y, z, t) can, by (2.2), 
also be regarded as a function of (x, y, p, 1) 
and we need a notation for the partial deri- 
vatives of fin the two cases. When fis regarded 
as a function of (x, y, z, f) we shall denote its 
We geen 

CL ON UE LS = 
when f is regarded as a function of (x, y, p, 
t), its partial derivatives* will be denoted by 
fo fur fo» fi, Consider now a set of five dif- 
ferentials dx, dy, dz, dt, dp related by 


partial derivatives by 


pee een ss 


op 7 
Op = ay Ox + SE dt 


op 
dy 


* In SI, these partial derivatives are written 


of 7 of (2 | 
Fe 2 e p op’ \ oe |p 


Io GE IMeW Ir We a 


and by 
Oz = 2, 0x + 2, Oy + Z,0p + 2 Ot. 


Substitution of either of these expressions 
into the other must lead to an identity in the 
differentials and we thus obtain 


dp , op Ca Or PE 
ni en dy No aloe es 
2.3 
op LP is Pos wg 
; Dr An ‘ 


in which, of course, p is regarded as a function 
of (x, y, 2, f) and z as a function of (x, y, p, t). 
As to the operator for differentiation following 
the motion, we write, when f (x, y, 2, #), 


fof, 
dt — ot De 


(2.4) 


in which u, v, w are also regarded as functions 
Oey eo t):.büb wen f (is 9, put). we 
write 


Df 


Di = +tuk tv + of, 


(2.5) 


where, u, v, & are functions of (x, y, p, f) and 


DDP 
Eid (2.6) 


[0) 


In Sutcliffe’s terminology, a weather situation 
in which @ is different from zero, is said to be 
“developing”, whilst if & is zero, there is 
no development. 
If we identify f with p and use (2.3), it is 
easy to show that 
d 


—— =, 


dt 


so that @ is the result of transforming 


(2.7) 


dp (x, y, 2, t) 
dt 
by. 2.2): 

In the following pages we shall be con- 
cerned with the (x, y, p) system in which w, 
strictly speaking, plays no part being replaced 
by Dz/Dt. Nevertheless is it convenient to 
retain w and to write 


w(x,y,p t) = A+n2,4+V2,+ 02, (2.8) 


3. The kinetic metric and the fundamental 
equations 
In order to translate the general formulae 


of MI and MII to our special case we shall 
bear in mind the following identifications 


(3-1) 


RAL) ee x ee Sy, x3 = p, | 


©, J 


K=HUl =u, P=? = v, p=u3 


in accordance with the usual convention of the 
tensor calculus. If then K is the kinetic energy 
of unit mass of air we have (MI (4.4) ), in 
the (x, y, z) system, 

Ke 


{x24 y2 +224 [k (z+ a) 


+ Ixy—kxz}, 


ly] à + 
(32) 


where a is the Earth’s radius. In this expression 
for the kinetic energy the “‘centrifugal” part 
— proportional to w? — is neglected for it 
can be regarded, in the equations of motion, 
as being incorporated in g. Using (2.2), (2.8) 
we have, in the (x, y, p) system 


NH 


= = {98+ v2+ (+ uz Vey + OZ)? + 
+ [k(z +a) —ly] u + lxv — 
—kx (z,+ uz + v2, + ©2,) }. 


(3.3) 


By (3.1) the corresponding kinetic metric is, 
(MII par. 2), 


ee (: =: 2) (das) — 


—- S Gardidé + 2y_,dxtde), (3.4) 

qr=i 

where 

Vaa = 22 —kxz, 

vu = = (2212. +k(z+a)—Iy—kxz,}, | 
ae 

ya= (222 + Ix — kxz, }, \ (3.5) 

MN = (22:2, — EX 2s Ms 

Y= IH, Yrr=I + 2,3, 33 = 2p, 

YI2 = Lx lys YZ = Sy Spy Y23 = Sy Sp. 


NY 


VPN SEPT PON KS 
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Since 712, 713, Y23 are not zero, the system (x, 
y, p) is not an orthogonal coordinate-system, 
but is “skew”. We shall require the determi- 
nant, A, of the y, (q, r=1, 2, 3) whose value 


by (3.5) is 


It is essential that /\ + 0, a condition satis- 
fied in the atmosphere where z, cannot be 
zero. 


(i) We can now write down the expressions 
for the divergence and the equation of conti- 
nuity which, by (3.1), (3.6) above and 
MII (3.4) and (5.12), are respectively 


te ee ~ 
=< = a - Ux. Vy == On, (3.7) 
1D (oz,) ~ 
02, = +, + Vy + Gp =O (3.8) 


en 
> 
. 

SS 


Since the coordinate-system is skew the 
simplest way of finding the equations of 
motion is to use the Lagrangian method 
(McCoNNELL 1931). These equations in 
terms of the notation (3.1) are 


D {9K\ 9K __1 ap 
Dt (Fa) le a u, 
(q= 1, 2,3), (3.9) 


where F, is the (generalised) force, apart 
from the pressure-gradient. In the (x, y, p) 
system we have, obviously, 


MEIN NIE (3.10) 


whilst if gravity is the only external force 
acting, F, is derivable from the potential 
gz so that 


Fı= — 92, Fr=—gz, F3= —g2p. 


(3.11) 


Using (3.3) in (3.9) and also taking account 
of (3.5), (3.10) and (3.11), the equations 
become, after some calculation, 


Du Dv Do | 
1 — ER us 
Dey LENS TAN CRT 
— (g + b—ku)z,— kw + lv, 


Dv Do 
ery a 


—(g + b— ku)z, — In, 


Du L Dv Dom 
Y13 Dt 23 Dt | Y33 Dt oar 


iD 
Yız Dr + V22 


I 
= (& + gz) ib kuls, | 


where 
D2 Dz Dz Dz 
b = u: x. ES ~ ee 
Di Hape UBER, 
(3-13) 


We can solve equations (3.12), with the 
help of (3.6), for the relative accelerations 
and thus obtain the final form of the 
equations of motion in the (x, y, p) system, 
Viz: == 


ae a 
Dt 
NE 
+ = (> + gz), A| 
2» \o p (3-14) 
Dv MET 
De u=—g2,+ = (2 ga). (3-15) 
Do 2% 
Dt + 2, (lv — kw) - 
2 I 
= lu E en N) — 
I I 
a (12e) (2 + = 
2, ( „) 0 8 =p 
(3.16) 


where w is given by (2.8). 


(iti) The rate of change of vorticity following the 


motion may be found thus: — By MII 
(3.13) and (3.5) above, the components of 
(absolute) vorticity are 


12 Gy € Mew ET PTE = 


Ge — Lu + (wa), — (w2)y}, (3-17) 
p 
G=k+— 1m — (wep) + (w2x)p}, (3-18) 
p 
i = - {] + v,—u, + (w2,). — (W2x)y— 
<p 


— kz, }, (3.19) 


the indices 1, 2, 3 referring to the x, y 
and p components respectively. The rate 


of change of £; is then by MII (6.8), 
Dé; 


oye UNE sO Cy Ogee 
5 X; + Y; 
+ a O° = Sa 
az 
where by MII (6.7) and (3.10), (3. 11) 
above, 


Ni (he (Fa), 3 re) = 9; 


ter 
3 0); À Ed . 


Hence by (3.7) 


Y 


Dé; 3 De; 
De = a C3 (the + vy) + 


+ Cr Wx + TOME (3.20) 


Now in (3.19) write 


Gem —,d = (wz,).— (W2x)y — key, 


(3.21) 


and then (3.20) is 


/ I <r IL. N 
an: Hd) ee = 
__ 1 Dep 5 =. 
Eon (UFER) 


or, finally, by (3.17), (3.18) the rate of 


change of the quantity (1 +¢+-¢') is given 
by the vorticity-equation 


IHN) = Un tm) (LE +0) — 


— 0,0, + Op + Oy (Wy Z_—Wa Zp + Re) — 


(3.22) 


— Wx (Wy Zy—Wy Zp) 
where (u,+ vy) is the isobaric divergence. 


(iv) Lastly the equation of heat transfer by bodily 
motion (conductivity, viscosity, radiation, 
etc., not explicitly analysed) is, by MII 
(7.5), 


DRAP ND DER 
Jp dés ae yi (=), (3.23) 


where J is the mechanical equivalent of 
heat, q is the heat content of unit mass of 
air and c, the specific heat of air at con- 
stant volume. 


It should be noticed that the equations found 
above are quite general within the limitation 
imposed by the use of local rectangular coor- 
dinates. In particular, they do not depend on 
the assumption that the equilibrium of the 
atmosphere in the vertical is hydrostatic. 


4. Dimensionless variables. Hydrostatic 
approximation 


We proceed to find approximate forms for 
the fundamental equations of the preceding 
section by an appeal to observation. To do 
this we introduce dimensionless variables, 
indicated by capital letters except in the case of 
the time, the density and the temperature when 
lower-case Greek letters will be used. We 
write, for the dimensionless variables, 


ro To zo 
a 
p CES à Le 
P=-,0=-,d=—. 
Po 00 31 


where po, 00, To are the values of the pressure, 
density and temperature at the point whose 
coordinates in the (x, y, 2) system are (ro, ro, Zo). 
We also write 
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DX mi V — Pb ait 
Dt se Se De? | fre” 
DP © D 7 w 
= He W= — = — 
Dr Ipo’ Dr lzo 


The conversion of some of the principal func- 
tions to the dimensionless variables is given in 
the following table. — 


no = LR = Zy; | 
Yo To 
Du _ „DU Dv _„,, Dk 
Dt cad by DR ADT 
Do Sig DIT. 
Dire Dr: 
DZ, 
b= beB= be. (= 
& u 
DZ .DZy DZp 
a U | / | . 
Dee Dr = De): 


, lzo2 
= = | (W2y)x— (WZx)v} — 
ee aay ae 
ro 


> lzo2 
"à pero 

(w 2x) — (w2,)x = 
D 1202 


Poto 


| (WZy)p — (WZp)y } ; 


{ (WZx)p — (WZp)x}. 


It is a fact of observation that the isobaric 
surfaces in the atmosphere are only slightly 
inclined to the horizontal, a gradient of 200 
ft. in so miles being large and 200 ft. in 100 
miles being appreciable. Thus we may say 
that z,, z, have values of about 7.6 x 1074 to 
3.8 x 10-4. If then we take z,, z, to be equal to 
5x 10-4 and assume that the corresponding 
values of Zx, Zy are of order unity, we have 


by (4.3). 


Zo 

© = ar , 
eat at (4.4) 
In temperate latitudes (say 9=60°) we also 
have 


g 
2 == (LOK (4.5) 


Consider the equation of motion (3.16) con- 
verted to dimensionless units by (4.3) so that 
it reads, after multiplication by 20/g, 


Zo l? BDIT ' Zx k zo 
g | Dr | Zp (v i To w)— 
Zy U | B À Zo? ZX + Ly 
Zp Zp To? Zp 
klro Ü I { 2 2 
2 Ni 2 (ae + Zÿ)| 
|» 


In order to find zo and ro, let us assume that 


Lea NE 1077 (4.7) 
g Fa 102 2 -5 > 4.7 
whence 
Parce kiro ke Zo | 
OR ern (4.8) 
Zola OC), fo 3, x LOW CI | 


The location of the “standard” point (ro, ro, Zo) 
is therefore at some 300 km from O and near 
the plane Ox y. Thus this point is not far from 
the stratosphere and the temperature there 
(To) is therefore low. However it is not 
necessary to know To at all accurately: its 
value must certainly lie between 200°A and 
300° A. In this range the coefficient of 1/o in 
(4.6) is, using (4.8) also, 


Po > RTo = 
£OoZ0 82 2a) 


(4-9) 


Thus, if in (4.6) we neglect all terms whose 
coefficients are of order 1074 or 1077 and re- 
tain those with coefficients greater than or 
equal to unity, we find 


14 GC. MevITTTEe 


rise u Zp = 0, 
20020 © 
or 
I 
a + g2, = 0. (4.10) 


Reference to (2.3) shows that this equation in 
the (x, y, z) system is 


Op FE 
ae I 


the familiar hydrostatic equation. Thus the 
hydrostatic approximation may be deduced 
from the observed fact that the isobaric sur- 
faces in the atmosphere are, very roughly, 
parallel to the Earth’s surface. 

The equation of motion (3.14) may be 


written, using (4.3), (4.7) and (4.10), 


Neglecting the term in this equation whose 
coefficient is of order 1074, the equations of 
motion (3.14), (3.15) may be written in dimen- 
sionless units as 


V = Zx + (4.11) 
DV 
Oe = Saree (4.12) 


The equation of continuity (3.8) reduces, by (4.3) 
and (4.10), to 

Og Vy =p = 0; (4.13) 
a form due to Godart (quoted in SI p. 371). 
The vorticity equation (3.22) becomes, using 


(4.3) and again neglecting all terms with 
coefficients of orders 10-4 and 1077, 


5 (t+ 8) =—(Ux + Va) (1 + 2) — 


—TIlx Vp + Ily Up. (4.14) 


The equation of heat-transfer we shall deal with 
in par. 6. 


Equations (4.11) to (4.14) are therefore the 
simplified forms of the fundamental equations 
which result from an approximation equiva- 
lent to the hydrostatic one. In the next section 
we shall consider further simplifications arising 
from certain definitions but we note one such 
definition now. The geostrophic wind is defined 
to be a wind whose relative acceleration 
(DU/Dt, DV/Dr) is zero. Hence by (4.11), 
(4.12) 


U =—Zy, V= 2x. 


If the wind is not geostrophic, the so-called 
“ageostrophic components” are, by (4.11), 
(4.12), (— DV/Dr, DU/Dr). Thus, apart 
from sign, the “‘ageostrophic wind” is merely 
another name for the relative acceleration of 
the true wind. For brevity we write (4.12), 


(4.11) as 


Usz— Zy—V, (4.15) 


V = Zx +UÜ. (4.16) 


5. The vorticity equation 


In the previous section we have seen how 
we can reject certain terms in the fundamental 
equations by means of an empirical approxi- 
mation. In the course of doing this we have 
also cleared the equations of the parameters 


g, I, ro, Zo, po, Qo and To so that we are now 


no longer able to perform further approxima- 
tions of this kind. We are therefore thrown 
back onto definitions of certain types of mo- 
tion which may, or may not, be realised in 
nature. An analogy may make this point 
clearer: in elementary hydrodynamics we 
begin with the study of the irrotational motion 
of an incompressible fluid, i.e. motions in 
which, by definition, a velocity potential ® 
exists such that 

ID 


ee es 


ID 


D 
— we —, 
dy dz 
But in writing down these three “restrictive 
conditions”, as we shall call them, we do not 
imply either that all fluid motions are irrota- 
tional or even that they are approximately 
irrotational: we are merely defining the 
simplest case that arises. 
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Turning now to the vorticity equation 
(4.14), the first restrictive condition we con- 
sider is expressed by the equation 

ITx Vp — 11y Up = 0. (5.1) 
It would be satisfied if one or other of the 
three following situations arose in an atmos- 
pheric motion: — 


(a) Up, Vp were negligibly small or zero, i.e. 
if there were no change of wind-speed 
or direction on passing from one isobaric 
surface to another, keeping (X, Y) fixed; 


(b) 17x, 11; were negligibly small or zero, ice. 
if there were no change of JJ on passing 
from one point to another in an isobaric 
surface; 


(c) Ix, Ty, Up, Vp were large but happened 
to satisfy (5.1) throughout the region of 
the atmosphere considered. 


Whether plausible or not, (5.1) is usually 
asserted in discussing the vorticity equation 
which then becomes, using (4.13) also, 


D. (+8) = — (Ux + Vy) (1 + 2) = 
= [Tp (t + 5) (5.2) 


The second restrictive condition we con- 


sider concerns the “ageostrophic” wind (—V, 


U). We have two mutually exclusive possibil- 
ities, viz: — 


(A) 


The ageostrophic wind is comparable 
in magnitude to Zx, Zy; 


The ageostrophic wind is negligibly small 
compared to Zx, Zy. 


(i) Consider first case (A) in which, by (4.15), 
(4-16), (4-13), 
Vx—Uy=Zxx+Zyy+ Uxt+ ia 
=V?Z+Ux+Vy, | 
Hp = — (Ux + Vy) = Vx— Uy. 


lun] 
In) 
dd 


I 


| (5-3) 


Hence, in this case, Æ + V:Z and 
IIp + 0, and we may write (5.2) as 


A 


oe (log (nr + 2) } 


= IIp. 


or, 


15 


This equation may be integrated with 
respect to P from an isobaric surface 
where P = Pı to another surface where 
P=P, (< Pi), and the result divided by 
P;—P;. We obtain 


P: 
I D ER 
poof (es PE) jap 
Py 
IN, D 


If we define a pressure-average, &, of 
E by 


P: 
I D 
tar 
Pr 
= pm, (leg (1 + 4) }, 


we obtain from (5.4) 
D | 
Dr 


returning to our original variables, and 
writing € = 18, 


| 
Dt 
which is a form of the Rossby potential- 


vorticity equation for the case of ageo- 
strophic motion. 


a= pao 


1+¢ 
Pas: 


(5.5) 


Consider next case (B) in which wo; 
V-= 0, $0that, “by=(41s), (4.70), (5-3), 
U=—Zy,V= AN ENT 202) lip — 0: 


(5.6) 


Using these expressions in (5.2), we have 


Lie ON Ee BW A 
== LENZ eZ pr (5.7) 


From this equation we can obtain a par- 
tial differential equation of the third 
order in Z by satisfying (5.1) in either 
of the following two ways: — 
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(1) Up=0, Vp=0, which by (5.6) lead 


to Zyp =0, Zxp = 0, and so to 
V7? Zp =0. This restrictive condi- 
tion is therefore equivalent to satis- 
fying (5.1) by means of (a); 


(2) IT= 0, so that ITx=0, Ily=0 and 
(5.1) is satisfied in virtue of (b). 


In either case (5.7) is 

V?Z— Zu N Zx + ZxV: Ay = © 
(5-8) 

which may also be written by (5.6) as 


V (1 + 


V2 Z + U(r + E)x 4 =, 


(5.8a) 


If we convert this equation to our original 
variables, use (4.7) to eliminate ro?, and 
write h = gz, we obtain 


I 


u (1+ C)y + v (1+ 6), + 


(hie ac hyy) = 0 
(5-9) 


which, in our notation, is Sutcliffe’s 
equation SI (11) for the case he describes 
as being that of ‘‘non-development” 
dp/dt =o (i.e. I1—0). The equation how- 
ever is valid also under the conditions 
Up = 0, Vp = 0, IIp = 0, which imply that 
U, V, II are functions of (X, Y, t) but not 
of P. Equation (5.9) can thus be estab- 
lished for a geostrophic wind provided 
that (5.1) is satisfied in virtue of (r) and 
(2) above. 


To obtain the remaining formulae of 
Sutcliffe’s theory we return to case (A) 
in which, by (5.3), & + V2 Z and 
IIp=# 0. Equation (5.2) may be written as 


V2 Ze + Uxe + Vyr + UEx + VEy + 
+ II&p = Ip + IDE. (5.10) 


If we now introduce the restrictive condi- 
tion expressed by the equation 


Uxr a5 Vie = ITp& — IT2p, 


(5.11) 


we find that (5.10) becomes, using (4.13) 
again, 


V2 Zp + Ux + VEy = lp = — 
— (Ux + Vy). (5.12) 


On returning to our original variables, 
with h= gz, we have 


hy) —{u(l +E) + v (ltd), = 


= Lu, + vy), 


us 
(5.13) 


which, in our notation, is Sutcliffe’s 


equation SI (10) for the case of ageo- 
strophic motion. It is therefore valid 
whenever it is known a priori that (5.11) 
is satisfied in the atmospheric motion 
considered. 


Sutcliffe however proceeds further by means 
of a second restrictive condition. Consider the 
equation (5.12) at two isobaric surfaces P—P: 
and P=P; so that 


Hence by subtraction 
(VV? Zz)2— (22): + UE x)2 — Ux(Ex)s + 
A ON wai Valeo ls eee 
(5.14) 


If we now impose the restrictive condition 
expressed by the equation 


(V? Z)1+ Ua(E x): + Va(By)1 = (V2Z,)a + 
+ Ui(E x): + Vi (Ey)», (5.15) 
and use it to eliminate (VW? 2): — (V2 Zi 


from (5.14), we can write the result in the 
orm 


(U,— Un) (1 = Ex SP 22) x ar 
ci (Vi — V1) (1 + 814+ 8,)y= (IIp)2— (ITp)1 = 
Te (Ux a Vy). Se (Ux + Vy): (5.16) 


Returning to our original variables, this 
equation is 


Ga 


5 


- 
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(12 — tx) ( m Cx = Ca)e — 
(nn) 4040), = 
I { (ug + Vy)2— (ux + v,)ı}, (5.17) 


which is equivalent to Sutcliffe’s final equation, 
SI (14), as may be seen by writing 


V = (12 — U1) i SE (v.—vı) j, ae ox,0= Cs, 


(WV. 7p) (J+ 64+ Lo) = (m—un) (1+01+02), + 
ed (et: + Sa), 


div, V = (u, + vy)2, div, Vo = (u, + vy). 


Thus equation (5.17) is valid only if it is 

own a priori that both the restrictive condi- 
tions (5.11) and (5.15) are satisfied in the at- 
mospheric motion under investigation. In 
terms of our original variables, with h = gz, 
these conditions are 


D D 2 = 
(5). Li ee =1(6,0—@6,), (5.114) 
E (han + Ayla + ua Cals + ve = 


= fh + Ada + ts (Can + 9% (Co 
(5.15 a) 


Hence Sutcliffe’s theory imposes more condi- 
tions on the motion than that which led to 
the potential-vorticity equation (5.5). It is 
true that each side of the equation (5.15 a) 
bears a superficial resemblance to the left- 
hand side of (5.9), which is valid for geostrophic 
motion. But Sutcliffe is not considering 
geostrophic motion, for the purpose of (5.17) 
is to deal with cases when (u, + vy), i.e. his 
div, V, is not zero. It also does not seem possible 
to deduce one of the equations (5.11a), (5.15a) 
from the other. Before, therefore, any synoptic 
conclusions are drawn from (5.13) or (5.17) 
it would seem necessary to investigate empir- 
ically to what extent, and in what synoptic 
situations, (s.11a) and (5.15a) hold good. 


6. The equation of heat-transfer 


The continued recurrence of \7?Z, in the 
formulae of the preceding section suggests 
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that further light would be thrown on the 
problem if an independent method for cal- 
culating Z, as a function of (X, Y, P,t) could 
be found. To take a simple example: in the 
geostrophic case, Z is a solution of the partial 
differential equation (5.8) and Z would there- 
fore involve a number of arbitrary “functions 
of integration”. Hence an independently calcu- 
lated value of Z, would serve as a kind of 
“boundary condition” and fix some of these 
unknown functions. Sutcliffe and Forsdyke in 
SII therefore seek such a determination of 
Z, by means of the hydrostatic equation (4.10), 
the gas equation p= Ro T and the equation 
2 Baer (3.23). We have in fact by 
4.10 


(6.1) 


epee ks 


and we also have the identity © = p,;+upy+ 
+ vpy +p, which, by (3.10), gives p, =o. 
Thus differentiating (6.1) partially with res- 
pect to f, we obtain 


LZn =—R 25 


5 (6.2) 


Integrating with respect to p from the isobaric 
surface p= pr to the surface p= p2, we have 


pı 
hy = ¢ | (z)2—(z)2} = RS Tid log p. 
pa 
(6.3) 


The equation of heat-transfer (3.23) and the 
gas-equation yield 


DT 
J& Dt 


©, 


Dq 
ID: 


whence, if J’ is the adiabatic lapse-rate, we 
obtain 


ie Er el 
T, = — (uT,+ vT,) + à ( r,) 
La 
py Lt 


If this is introduced into (6.3) we have 


18 Gy € MeV = 


5): 
(6.4) 


which is, in our notation, the equation SII (3). 
In order to find the magnitude of its various 
terms, we convert to dimensionless units by 
(4,1), (4.3). and introduce also the dimen- 
sionless heat-unit Q = q/qo. We find 


Py 


RATIO 
Aen f{-Wox+ Voy) + 
4 oy, 
Ber qo = 
| = [ex 
H(t 5 0) } eatin De d log P. 
(6.5) 
Since 
UR 2508 Mo A 
© 6.2405, 2 4.1863 X 101 “10° 


it follows that all the terms in the integrand 
of (6.5) are of the same order of magnitude 
provided that the heat-content, go, and the 
temperature, To, of unit mass of air at our 
“standard” point (ro, ro, Zo) satisfy the condi- 
tion go = € To. If therefore the right-hand 
side of (6.5) can be evaluated directly, we shall 
have a value of Z, at an arbitrary level for 
use in the equations of par. 5. 

We note however that a study of equation 
(6.4) by means of the data of observation 
does not, by itself, throw light on the dynamics 
of an atmospheric motion. For this equation 
depends only on the hydrostatic equation and 
the equation of heat-transfer. It is only when 
(6.4) is used in conjunction with equations 
derived from the two “horizontal” equations 
of motion and the equation of continuity that 
dynamical consistency can be hoped for. 


7. The equivalent-barotropic atmosphere 


The methods of this paper can be applied 
to the analysis of a form of the vorticity-equa- 
tion which occurs in the work of J. G. Charney 


(1949). This author employs the (x, y, 2) 
system in his work but, in order to permit com- 
parison with our results, we shall convert his 
formulae to the (x, y, p) system. We deal 
again with case (A) of par. 5 and once more 
assert (5.1). Following Charney, the first 
restrictive condition on the motion is expres- 


sed by 
U =A(P)U'(X,Y,7), V=A(P)V' (XY, 0), 


(7-1) 


where A is an unknown function of P. This 
condition means that, for fixed (X, Y, 7), 
there is no change of horizontal wind-direc- 
tion, but a possible change of wind-speed, 
as we pass from one isobaric surface to another. 
Let P =P: now correspond to the Earth’s 
surface whilst P = P; = o corresponds to the 
“top” of the atmosphere where we also 
assume that JJ = 0. We define P= P to be an 
isobaric surface on which A(P) is equal to its 
pressure average, i.e. 


Let 


so that, if U, V refer to some isobaric surface 


other than ee 12% we have 


À — À = 
Uso ver NOR 
4 | 
We define 
ia = Vz oe U; > 
so that 
See Vy Uy = = Say TA 
2 (4 
Hence the vorticity equation (5.2) becomes 
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since Zp = 0. We next define the dimensionless 
number 


K={ j'A(PdP\| [fA)dP], (76) 


9 


which, apart from a reversal of limits in the 
first integral that is probably due to a misprint 
in his paper, plays the same part as Charney’s 
K. If we then integrate (7.5) with respect to P 
from o to Pı and use (7.2) and (7.6), we find 


ergs 
E,+ K(USx + Vey) + = [ Hardp = 
ie} 


A 


We now impose a second restrictive condition 


expressed by 


(7) 


Pi 
[ ITA» dp= O, (7.8) 


which corresponds to Charney’s neglect of 
terms in the vertical velocity-component 
w. Considering also the equation of conti- 
nuity (4.13) we obtain, by (7.2), 


or, if IT = II: at the Earth’s surface 
He (Ux an Vy). 


Hence integrating (7.8) by parts, using the 
equation of continuity and this formula for 
IT, (+ 0), we find 


APE RAP) =», (7.9) 
which shows that the restrictive condition 
(7.8) is essentially a limitation on the func- 
tion A. Again 


ar 


Thus finally (7.7) may be written in either of 
the equivalent forms 


(7.10) 


Ale +K(UEx+ Vey) — IH(r + KE) = (©. 
(rt) 


Returning to our original variables by (4.3) we 
have 


poA (Pips) = Ÿ Alpip)dp 


pr 


Pi 
K={pof Ardp} |(S Adp)», 
and (7.10), (7.11) become, respectively, 


G4K @l, +96) ++ KO m. +2,) =0, 


(7.12) 
ogg (E+ KGET + TE) =o, 
(7.13) 
where 
a= IroU, 
5 = IrV, 
Die 


The equation (7.12) corresponds, in the (x,y,p) 
system and for / = constant, to an equation 
given by Charney (1949, foot of page 383). 
However, (7.13) differs from his next equa- 
tion (1949, top of page 384) in the following 
respects:—By (2.7) and the hydrostatic equa- 
tion we have, for the last term, 


20 GACAMEVIIANNTE = 


Or IS agp oP 
Spe a Laer oF 
an 
i ay Eee 


whilst Charney has — in our notation — 


Ox i op ne \ 
po un pi { Di gow = 


The discrepancy in sign is explained by an 
error due to a reversal of the limits of integra- 
tion in Charney’s manipulation of the tendency 
equation at the foot of page 383. The substitu- 
tion of the surface pressure pı for our po is a 
consequence of the use of local rectangular 
coordinates. The pressure po is that at our 
standard point (ro, ro, Zo) which is near the 
plane Oxy. In the “fat earth” theory which 
results from the use of local rectangular coor- 
dinates, po could be identified with pr, though 
we have preferred not to do so in par. 4. 
Apart from these two points, our formula 
agrees with Charney’s in the case of geo- 
op op 


strophic motion when # — 


ox dy 


We return to the consideration of the geo- 
strophic case, and write (5.8 a) as 


(1+ 2),+U(1+8)x+V(1+8)y =0 


’ 
or, returning to our original variables by (4.3), 
we have 


7+), +u(l +8), +v(1+ 0, = 0. 
(7.14) 


This equation corresponds, in the (xy ap) 
system, to the equation [1949, (42)] taken by 
Charney as governing the motion of the 
equivalent barotropic atmosphere. Since, by 
(7.1), Up, Vp are not zero, the condition (1) 
of par. 5 (ii) is not satisfied and therefore 
condition (2) must apply. This means that 
(7.14) is valid when &=o (Sutcliffe’s non- 
development) throughout the atmospheric 
motion, a restriction which replaces Charney’s 
condition that op/ot shall be zero at the 
Earth’s surface. Thus in the geostrophic case, 
Sutcliffe’s equation (5.9) is identical with 
Charney’s (7.14) and both are equivalent 
forms of the third-order partial differential 
equation (5.8) for Z. 

If the foregoing analysis is accepted, it 
appears that Rossby, Sutcliffe and Charney 
are each using models of atmospheric motions. 
All three models use the restrictive condition 
(5.1) and, in Rossby’s case, no further restric- 
tions are needed in order to obtain the poten- 
tial-vorticity equation (5.5). This appears to 
be the most general case for ageostrophic 
motion. At the other extreme is Charney’s 
equivalent barotropic atmosphere, with the 
additional restrictive conditions of geostrophic 
motion and (7.1), (7.9). Whether these addi- 
tional restrictions are accepted or not, they are 
certainly easy to interpret in physical terms. 
Intermediately comes Sutcliffe’s model, in 
which the additional restrictive conditions are 
expressed by equations (5.11 a) and (5.15a). 
Whilst Sutcliffe has thus escaped from the 
case of geostrophic motion it is not apparent 
that his additional restrictions have any ob- 
vious physical interpretation. 
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Abstract 


Some general comments are given regarding numerical methods of integrating the 
hydro-dynamic equations. A series of tendency computations for Europe and eastern At- 
lantic are presented, by and large following the scheme given by BoLIN and CHARNEY 
(1951). Data requirements for this kind of computations are discussed on the basis of the 
difference in the success of the method over regions with poor and with good data. 


During the last few years several groups of 
meteorological research workers have devoted 
an increasing amount of effort to the problem 
of translating basic dynamic principles into 
practically usable methods for quantitative 
objective forecasting of changes in the large- 
scale flow patterns of the atmosphere. The 
advent of high-speed electronic computing 
machines has been a factor of decisive impor- 
tance in this development; the fact that we 
now at long last seem to know how to construct 
theoretical models of the atmosphere which 
combine simplicity with the ability to repro- 
duce several highly characteristic features of 
the observed atmospheric motions has un- 
doubtedly contributed to the growing interest 
in the development of forecasting methods 
based on dynamical principles. 

The most systematic approach in this field has 
been made by the group working at the Insti- 
tute for Advanced Study in Princeton under 
the direction of J. Charney and J. von Neu- 
mann. They started with a very simple model 


1 The investigation described in this report has 
been supported primarily by a grant from the U.S. 
Weather Bureau. 


of the atmosphere by assuming autobarotropy 
and two-dimensional flow, but retaining the 
principle dynamic consequence of the spherical 
shape of the rotating atmospheric shell. With 
this model as a starting point they were then 
able to build up more complex structures capa- 
ble of accounting for atmospheric features not 
contained in this initial model. After prelimi- 
nary tests with a one-dimensional model (CHAR- 
NEY and ELIASSEN, 1949) the results of the first 
attempts to make complete 24-hour forecasts 
using the two-dimensional model were pub- 
lished little more than a year ago in an article by 
CHARNEY, FJÖRTOFT and von NEUMANN (1950). 
As one might expect the results were by no 
means perfect but at least very promising. 
However, only four complete forecasts were 
made; high speed electronic computing ma- 
chines and dependable soo mb charts, carefully 
analyzed with reference to the vorticity distri- 
bution, must be available before a systematic 
testing program can be carried out. 

If the hydrostatic and geostrophic assump- 
tions are introduced in the equations of motion 
as proposed by CHARNEY (1948) the system of 
differential equations may be reduced to one 
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equation which is only of the first order with 
respect to time, in other words the changes of 
the flow to be expected are completely determined 
by the state of motion of the atmosphere while no 
knowledge of the rate of change of this state is 
necessary. It is this fact that provides the basis 
for all the numerical integrations through 
iterative processes which have been proposed 
until now. Regardless of the ultimate prognostic 
value of this fundamental result it is evidently 
of great scientific interest to explore its validity 
further. It is in this connection of particular 
interest to know to what extent this is true for 
the still more simplified barotropic model of 
the atmosphere, e.g. to determine the extent 
to which the observed motions in the soo 
mb level, interpreted as if they were motions 
of an equivalent barotropic atmosphere, deter- 
mine the short range evolution of the circula- 
tion at that level. For this reason the Princeton 
group working under Charney’s direction in 
1951 initiated a series of tendency computa- 
tions from selected soo mb charts and the 
results of these were recently set forth in a 
paper by Borın and CHARNEY (1951). 

Before proceeding to a discussion of the 
results of the series of similar computations 
carried out here a few further remarks should 
be made. 

In any attempt to develop objective methods 
for forecasting {61S extremely important to 
keep in mind the limitations that are prescribed 
by the amount of observational data that is 
available as well as the accuracy of these data. 
It is already obvious from the previous investi- 
gations that not even the two-dimensional 
barotropic model in its present form can be 
used to its fullest extent over the oceans where 
the aerological network is sparse. It would 
then be very unrealistic to try to apply a de- 
tailed three dimensional model of the atmos- 
phere over such an area. On the other hand we 
know from experience that it is possible with 
existing subjective methods to forecast at least 
some large scale phenomena even with the 
small amount of data that now is at our disposal. 
As these large-scale processes in the atmosphere 
undoubtedly to a large extent are barotropic 
in their character it seems that one of the most 
important problems at present is to develop 
methods that make it possible also to forecast 
such developments objectively. In such a 
method it should then in principle not be nec- 


essary to consider the’ detailed structure of the 
flow as is done at present since such details are 
not known. 

There is another reason for studying the two- 
dimensional barotropic model in order to 
utilize its potentialities more fully before at- 
tempting to integrate the three-dimensional 
vorticity equation. Such an integration means 
such an extension of the computational pro- 
gram that there would be great difficulties to 
utilize the results in practice, due to the lim- 
ited capacity of present electronic computers. 
From a practical point of view it therefore 
seems as essential to improve this barotropic 
model as to attempt to introduce new phys- 
ical factors which, though important, would 
also increase the computational work con- 
siderably. 

As is seen from the basic equation for the 
barotropic model the computations involve 
the evaluation of vorticity advection, which 
means a differentiation of the initial contour 
field three times. It is obvious that no great 
accuracy can be obtained for such a quantity, 
in particular if no objective method of smooth- 
ing the original contour field is used. However, 
these inaccuracies are of relatively small con- 
sequence for the accuracy of the final tenden- 
cies because an integration is performed twice 
in order to obtain them, which means a 
smoothing and elimination of random errors. 
From a practical point of view then the ques- 
tion arises: Is the procedure of first differen- 
tiating three times and then integrating twice 
necessary in order to obtain the tendency field? 

In view of what has been said above and the 
benefits that would result from the successful 
development of an adequate technique for 
numerical prediction it was decided to make 
this one of the principle fields of research at the 
expanded and reorganized Institute of Meteoro- 
logy of the University of Stockholm. As a first 
step it was decided to perform a series of baro- 
tropic tendency computations of the type re- 
ferred to above, in order to enable the various 
members of the research group to become 
familiar with the technique of the computa- 
tion methods. Future reports will deal much 
more fully with these investigations. At the 
present time we shall restrict ourselves to a 
brief preliminary discussion of these tendency 
calculations, which for the first time provide 
us with some information concerning the 
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applicability of the barotropic model in North- 
western Europe. 

The active work of our group is at present 
supported by grants from the Swedish Natural 
Science Research Council, the Knut and Alice 
Wallenberg Foundation, the University of 
Chicago and the U.S. Weather Bureau. We are 
also deeply indebted to the various organiza- 
tions in Sweden and abroad who have contri- 
buted to our project by making professional 
personnel available. The participants in the 
tendency computations reported here have 
been: Mr G. Arnason, Iceland; Mr B. Bolin?, 
Sweden; Mr Ph. Clapp, USA; Dr A. Eliassen}, 
Norway; Mr K. Hinkelmann, Germany; Mr 
E. Hovmöller, Denmark; Lt. W. Hubert, 
USA; Dr E. Kleinschmidt, Jr., Germany; Dr 
C. W. Newton!, USA; Mrs H. Newton, USA; 
Dr H. Schweitzer, Germany; Miss Ch. Steyer, 
Germany. 


The method used in these computations is 
the same as that used at Princeton and described 
by Bouin and CHARNEY (1951). Computations 
were made on a rectangular grid, with grid 
interval of approximately 310 km at latitude 
50° N. Because of time limitations, computa- 
tions of vorticity, vorticity advection and ten- 
dency were made only at alternate points on 
the grid, in the manner discussed by Bolin and 
Charney. It was necessary to restrict the area 
over which the computations were made to a re- 
gion about 4,000 km square, with its center near 
the British Isles (see figure 6). An important 
part of this region lay over the Atlantic Ocean 
so that errors both in computation and in veri- 
fication arise from the uncertainty of analysis 

‚over this region with sparse observations. 

Two methods of verification of the com- 
puted tendencies were used. The first of these 
is the method used at Princeton, that is, to 
compare the computed value of 


I Oz Oz 
2 Be R aa 


with the observed height changes 


Xo = Zt=hb+A — “t= too 


1 Project leaders. 
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At being 12 hours. z is the height of the soo mb 
surface. In all, 14 pairs of tendency computa- 
tions at the soo mb level have been made here. 
As pointed out by Bolin and Charney, verifica- 
tion by the above method tends, for fast- 
moving systems of period less than 36 to 48 
hours, to give poor results which are not related 
to the efficiency of the computation method 
itself. In the second method, use was made 
of the frequent upper-air observations from 
Great Britain, in a manner to be discussed 
later. 

Since in solving for the tendencies it is 
necessary to make some assumption regarding 
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Fig. Ib 


Contours (solid lines) and computed tendencies (dashed lines), (a) 1500 GMT 18 November 1951; 


(b) 0300 GMT 19 November 1951. Contours labelled in tens of meters, interval 40 m, tendencies in tens of 
meters per 12 hours. In wind symbols, solid triangles represent so knots, full barbs 10 knots. 


the tendencies at the boundaries, the verifica- 
tion procedure was applied only in the central 
part of the region, excluding that part (two 
grid points from the boundary of the region 
over which computations were made) which 
is strongly affected by boundary conditions. 
Computations were made first on the assump- 
tion that the tendencies at the boundary are 
zero, and then repeated using, for each pair of 
computations, the actual 12-hour changes at 
the boundaries. The results of the verification 
are presented in table I. x, denotes the value 
of the mean tendency when it is assumed that 
the tendency on the boundary is zero and x, 
is the value when the observed values x, are 
used as boundary values. A bar indicates a mean 
value over the area of verification. r, is the 
correlation coefficient between x, and xy, ra the 
correlation coefficient between x, and x. 


a = V2' (x, — x)2/n, ie. & gives the root- 
mean square of the error in x, and correspond- 
ingly &, gives the error of x9. Finally o, denotes 
the root mean square of the observed changes 
(x9), in other words o, is the error of a forecast 
of no change. 

The following should be noticed: 

1. From all cases combined, an overall corre- 
lation coefficient of 0.7 was obtained, somewhat 
lower than that reported by Bolin and Charney. 
The difference is hardly significant. On the 
other hand these values on the average are con- 


siderably higher than those given by SAWYER 
and BusxBy (1951) for two computations. It 
seems safe to conclude from the computations 
made at Princeton and those reported here that 
the results of Sawyer and Bushby are not re- 
presentative for computations of this kind. As 
seen from the table poor correlation was ob- 
tained in two cases. In both these cases (as well 
as in some of the others with fairly low correla- 
tion) the continuity of the analyses was not the 
very best. To a certain extent this depended 
upon the fact that the analyses were made by 
different members of the group. However, the 
last ten maps (Dec. 5, 15—Dec. 10, 03) were 
analyzed by one of us with special attention 
to continuity in particular over the Atlantic 
Ocean, where the data are not sufficient to 
determine the flow patterns completely. Here 
the correlation is considerably higher. It should 
be noted that in no case correlation below 0.5 
was obtained in the computations at Princeton, 
which also suggests that such low values are 
mainly due to inadequacy of the initial data. 
In those cases all computations were made over 
the North American Continent where quite 
a dense net of radiosonde observations is 
available. 


2. By and large there is a very little difference 
between the values obtained by using x, or x, 
in the verification, which indicates that the 
boundary influence in most cases is negligible 
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AVERAGE COMPUTED 
TENDENCY 


IBISZ -1911037Z 1951 


OBSERVED CHANGE 
18-11-15Z - I9-11-03Z 1951 


Fig. 2b 


(a) Average computed tendencies; (b) observed 12-hour change, 1500 GMT 18 November 1951 to 0300 


GMT 19 November 1951. Values given in tens of meters per 12 hours. 


at a distance of more than two grid points from 


the boundary. 


3. In one case out of these fourteen computa- 


tions the root mean square of the errors was 
greater than that of the actual changes. It is 
quite obvious that the root mean square of the 
errors becomes quite large if the computed 
changes are systematically overestimated as is 
the case in these computations. Such an over- 
estimate is not necessarily a defect of the baro- 
tropic model but may depend upon our way 
of applying the theory to the real atmosphere 

y considering the flow at soo mb, or upon 
the geostrophic assumption. Furthermore it 
should be stressed that if the root mean 
square error is as large as that which would 
be obtained by assuming ‘no change’, this 
does not mean that the forecasts are of no 
value. If the correlation differs from zero in 
a significant way, the computed pattern still 
to a certain extent resembles the actual change 
pattern. Both the correlation coefficient and 
the root mean square errors are strongly af- 
fected if there are a few large errors. Thus if 
the result is good over 75 % or more of the 
area, with large errors in the remaining part, 
the root mean square errors will be quite large. 

4. In columns 6—8 of table 1 the computed 
and actual mean changes are given for the area 
of verification. There is by no means a perfect 
agreement but a definite correlation exists be- 
tween the computed and observed values. This 


is of considerable importance for determining 
large-scale changes of the flow pattern, quite 
separately from the changes associated with the 
smaller but fast-moving systems. 

Figs. 1— 5 show the results from two pairs 
of computations. The first one (Nov. 18—19) 
represents a case that is somewhat poorer than 
the average (r = 0.52) while the other (Nov. 
25—26) is one of the more successful computa- 
tions (r = 0.74). 

In the computed patterns Nov. 18—19 
(figs. 1a—1b) there is only fair continuity 
between the two tendency computations. 
The major failures in that respect usually 
depend upon inconsistencies between con- 
secutive analyses. For instance, the trough 
line across the North Sea November 19 
gives rise in the computations to a fall area 
to the north and a rise area to the south 
which did not exist twelve hours before. 
Therefore, a considerable difference exists 
between the observed and computed pat- 
terns in that area (fig. 2). Obviously we do 
not know the exact character of the flow 
as no observations are available from the 
North Sea, and the character of the flow has 
to be inferred from surrounding observations. 
This clearly shows the necessity of a dense 
net of homogeneous observations for present 
methods with a small grid size. 

With respect to the second pair of com- 
putations (figs. 3— 5) the following should 
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Fig. 3. (a) Contours (solid lines) and geostrophic absolute vorticity (short dashed lines); (b) contours and computed 
instantaneous tendencies (dashed lines), 1500 GMT 25 November 1951. Units as in fig. 1, vorticity in unit of 107° sec!. 
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Fig. 4. (a) Contours and geostrophic absolute vorticity; (b) contours and computed instantaneous tendencies, 
0300 GMT 26 November 1951. Notations as in Fig. 3. 
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OBSERVED ANGE 


Pe - 2611037 1951 


Fig. 5a 


Fig. 5. (a) Average computed 


Fig. 5b 


tendencies; (b) observed 12-hour change, 1500 GMT 25 November 1951 to 


0300 GMT 26 November 1951. Units are tens of meters per 12 hours. 


be noted: According to the analyses a pro- 
nounced increase of the absolute vorticity 
takes place to the northwest of Spain. This 


Big. 6. 
cients of observed vs. mean computed tendencies, for 
13 pairs of computattons. Observed geostrophic vor- 
ticities Were determined for each case utilizing heights 


Area of computation, and correlation coeffi- 


over entire region shown; instantaneous tendencies 
were computed for area inside dashed rectangle, and 
statistical verification was carried out for region inside 
smaller dash-dotted rectangle. Isolines of correlation 
coefficients are drawn for correlations made separately 
at 12 points inside smallest rectangle, each point in- 
corporating tendency comparisons at four surrounding 
points on computational grid, thus representing 52 
verifications of tendency pairs at individual points. 


result may, however, largely depend upon 
the analysis, the number of observations is 
insufficient to establish this indicated increase 
with any certainty. The agreement between 
computed and observed changes (fig. 5) is 
good, except in the northwestern part of 
the area, where the artificial boundary con- 
ditions influence the results considerably. 
Boundary assumptions are especially impor- 
tant in the case of large-scale systems, but 
often are significant also for the small-scale 
circulation patterns. 

In both of the examples given here (see 
figs. 2 and 5), the mean computed tenden- 
cies are for the most part considerably larger 
than the observed changes. This was in gen- 
eral true for the whole series of computa- 
tions, and most likely indicates that the 
soo mb surface is appreciably higher than 
the equivalent barotropic level. 

It was obvious at the outset that the tendency 
computations would be better over some re- 
gions than over others, merely from the fact 
that the density of the aerological network 
varies over the computational region. An at- 
tempt was made to obtain the geographical 
distribution of the correlation coefficient using 
thirteen pairs of computations (the one ex- 
cluded was done for a slightly different region). 
The results are shown in fig. 6. For each corre- 
lation coefficient the four surrounding points 
have been used, so that each of them was 
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derived from 52 computed and observed ten- 
dencies. Correlations are quite high near the 


‘center of the area (= 0.92); this is probably 


largely the effect of the reliable analysis due to 
the adequate coverage of good aerological sta- 
tions in the British Isles, but probably also 
depends upon the fact that this area is far away 
from the boundaries. The correlations drop over 
Oceanic areas, as would be expected, but also be- 
come significantly lower over western Europe. 
The inhomogeneity of the upper air material 
is undoubtedly responsible for at least part of 
this disagreement. 

One disadvantage in the method of verifica- 
tion used above is the fact that two computa- 
tions are involved in each verification. When 
computing the mean tendencies for the twelve 
hour period the random errors tend to decrease 
and thus the final results may be improved. A 
separate verification was made over the British 
Isles using the six-hourly observations that are 
available for eight stations in the area in order 
to ascertain that this fact did not play an im- 
portant role in improving the correlation coef 
ficients. No success was obtained in attempting 
to estimate the instantaneous tendencies at these 
stations from a time graph of the soo mb 
heights for the reason that these graphs showed 
short period oscillations of such a nature that 
the determinations could not be made objec- 
tively. Instead, the twelve-hour change in 
height at each station, using observations six 
hours before and six hours after the time of 
each individual tendency map, was compared 
with the tendencies read off at the proper loca- 
tions on the analyzed computed tendency 
charts. In this way it was possible to verify the 
computations for 28 separate map times. The 
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Toni. s Aa cn toners 0.89 99 70 
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results are shown in table 2, where r, is the 
correlation coefficient between computed and 
observed tendencies, 0, and o, denote the root 
mean square of these two quantities respec- 
tively. 

We notice that the correlations over this 
area average considerably higher than those 
computed for the area as a whole and they are 
in some cases surprisingly high remembering 
that only barotropic processes are considered. 
As seen from the second and the third column 
of table 2 the computed tendencies are about 
30% larger than the observed values. This 
may to a certain extent depend upon the fact 
that the observed values represent a mean 
tendency over twelve hours while the com- 
puted values are instantaneous, but this prob- 
ably does not account for the whole discrep- 
ancy. A certain overestimation of the ten- 
dencies seems to result from the assumptions 
used in the present method. 

Finally some further remarks might be made 
on the effect of analysis upon the tendency 
computations. The Princeton group found that 
little improvement could be obtained by re- 
analysing the series of s00 mb charts used in 
their computations. This was not surprising, 
since the area of their computations extended 
mainly over North America, where the ob- 
servations are quite dense and consequently 
large changes in the analysis could not be made. 
In the present series of computations, made 
largely over oceanic regions, it is evident that 
the analysis plays an extremely important role. 
For example, a re-analysis of the chart for 15 
November 0300 GMT, guided by the prin- 
ciple of eliminating large vorticity gradients 
introduced by the former analysis raised the 
correlation coefficient considerably for the 
corresponding pair of maps. It is quite probable 
that improvements in the verification could be 
obtained for several of the computations, if 
re-analyses were made. There is obviously a 
need for extensive study of the problem over 
regions with sparse data, particularly since the 
influence of the flow patterns over oceanic 
regions is rapidly and strongly felt over the 
neighboring land areas. In addition to the 
oceanic region, great difficulties are met in 
analysing charts over southern Europe and in 
the whole Mediterranean area, due partly to the 
fact that there has been no international stand- 
ardization of radiosonde instruments and 
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possibly also to errors in taking observations, 
as well as irregular reception of the observa- 
tions. Hence it is often quite difficult to arrive 
at a satisfactory picture of the flow pattern 
in this region. A definite need is also felt 
for an increase of radio wind observations, 
a field in which the British have an out- 
standing lead. 

In conclusion, we may let the data presented 
above and those given by Bouin and CHARNEY 
(1951) speak for themselves. Although nothing 
concrete can be said about the significance of 
the statistical quantities given above in regard 
to what will be obtained by further iterative 


forecasts made on high-speed computing 
machines, the agreement between computed 
and observed tendencies is certainly close 
enough to be very encouraging. This suggests 
that the effect of conservation of vorticity 
contributes for short periods of time by far the 
largest share to the observed pressure changes 
in the atmosphere. Other effects including 
baroclinic development of the type discussed 
by Surc irre (1947) are obviously of great im- 
portance for large-scale slow changes, and for 
more rapid changes in smaller regions of the 
atmosphere. 
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Abstract 


In treating a still limited material of observations on the content of nitrogen in pre- 
cipitation, collected under the supervision of Dr Egnér of the Agricultural College at 
Uppsala, we have obtained some preliminary results, which may be summarized as follows: 

The concentration of NH,-N and NO,-N in precipitation is dependent upon the 


amount of precipitation previously fallen. 


A mathematical expression is given for the dependence of concentration on amount 


of precipitation fallen. 


With aid of this relation a computation has been made of the 


normal supply of nitrogen at Uppsala, using known normal frequencies of different 


amounts of precipitation. 


A pronounced dependence of the concentration on the character of the air mass has 
been found to exist. Thus the tropical air contains between 10 and 30 per cent greater 
amount of nitrogen than the polar air, and almost double the amount found in the 


arctic air. 


From an agricultural point of view the 
question concerning the composition of pre- 
cipitation is an important one. On the basis 
of earlier available data it can easily be shown 
that the supply of plant nutrients in this way 
is not neglible. Thus the supply of nitrogen in 
the form of NH,-N and NO,-N in Sweden 
amounts to between 1 and 5 kg per ha!, a 
quantity much smaller but yet comparable 
to that generally furnished in an artificial 
way in the cultivated districts. 

The interest of agriculture in the question 
just mentioned was the chief reason why 
Messrs Torstenson and Egnér of the Agri- 
cultural College of Sweden took up the prob- 
lem in 1947 for closer examination. They 
applied modern methods partly devised and 
developed by Egnér. 

The problem has also, however, a geophys- 


1 ha is used as an abbreviation of “hektar”’ = 104 m?. 


ical side, which is not less interesting. The 
question regarding the presence and nature of 
condensation nuclei in the atmosphere is 
closely connected to the composition of the 
precipitation as regards dissolved salts. Their 
concentration and properties may be expected 
to throw some er on the difficult problem 
concerning the formation of clouds and rain 
in the atmosphere. Furthermore the composi- 
tion of the rain may show some interesting 
relations to the general properties of the air 
masses from which it is falling. Do the different 
air masses—tropical, polar or arctic air— 
show some caracteristic differences as regards 
the condensation products arising within them? 

On .account of the importance of these 
questions from a geophysical point of view, 
the present authors were attached to the en- 
terprise, at first to some extent for contributing 
to the organization of a synoptic net of meas- 
uring stations and later for treating the mate- 
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rial of observations with regard to meteoro- 
logical aspects. 

According to the program the technical and 
agricultural side of the matter would be treated 
by Dr Egnér and his collaborators. In their 
work would be included also the question of 
the annual and geographical variation of the 
composition of the precipitation. 

In the following we are adhering to this 
plan and will report on some theoretical and 
statistical results as regards the content of 
NH,-N and NH,-N in the precipitation in 
Sweden. 


I. The observations 


The samples of precipitation, considered in 
the following, were all taken at the Agricul- 
tural College at Ultuna (Lat N 59° 49’ Long 
E 17° 40’). The technical procedure is described 
in detail by Hans EGNÉR, ERIK ERIKSSON and 
ARNE EMANUELSSON (1949). Ion exchange re- 
sins were used for the analysis of NH,-N and 
NO,-N. With them the precipitation could 
be filtered practically immediately, thus lessen- 
ing very effectively the time of contact with 
dust, insects and pollen. Careful tests were 
made in order to clarify how completely and 
rapidly the ions were taken up by the filter. 
Special arrangements were devised in order to 
eliminate errors from contamination through 
bird dropping, pollen, insects etc. 

The observations cover the period 1. 10. 
1947—30. 9. 1950. The measurements of the 


composition refer to the different occasions of 
precipitation. Thus, values of the total con- 
tent of NH,-N and NO,-N in a given rain or 
snowfall of longer or shorter duration are 
obtained. The result of our analysis will show 
the desirability of a more detailed investigation 
of the changes in composition with time in the 
separate cases. of rain or snow. The observa- 
tions treated here refer, however, only to the 
integral content in the different cases. In the 
total, 117 cases of precipitation have been exa- 
mined. 


II. Analysis 


1. Change of concentration of NH,-N and 
NO,-N with duration of precipitation. 


Already, a superficial examination of the- 
values obtained show a rather marked depend- 
ence of concentration on duration. 

From Table 1 and 2, which show the fre- 
quency distribution of the concentration of 
NH,-N and NO,-N in milligrams per liter 
(mg/l) and its dependence on the amount of 
precipitation, it is evident that the mean con- 
centration of the total amount of precipitation 
is on the whole a decreasing function of the 
amount. It is consequently as if the precipita- 
tion drew its content of NH,-N and NO,-N 
from some source which is contributing less 
and less as the condensation and precipitation 
progresses. If we try to find a mathematical 
formula expressing the relation between the 


Table 1. Frequency distribution of the mean concentration of NH,-N (mg/l) for all air masses, 
given as a function of the total amount of precipitation 
Precipitation mm 
YNHi-N : 
dl 0.0] 2.0} 4.0] 6.0] 8.0] 10.0] 12.0] 14.0] 16.0] 18.0] 20.0] 22.0] 24.0] 26.0] 28.0] 30.0] 32.0] Z 


1.9} 3-9] 5-9) 7-9 
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0.00—0.19| 4 3 2 6 5 — 3 | — 2 I I | — | — I | — I 32 
0.20—0.39| 5 | Io 6 7 5 I 3 I IIl— | — I 40 
0.40—0.59| 4 Ole liées: 3 a I2 
0.60—0.79| 7 3 2 I 1 | - - 14 
0.80—0.99| 3 4 I | — I — | — 9 
1.00—1.19| I 3 — - 4 
1.20—1.39| 3 | — I | — 4 
1.40—1.59| — I | — | — a 
1.60—1.79| — | — I - I 
2} | 27 29 | 13) | 74 | 4] af 4] 2 I :|-|-| :|-| : 177 
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Frequency distribution of the mean concentration of NO,-N (mg/l) for all air masses 


mm 


Table 2. 
given as a function of the total amount of precipitation 
RReGi pil tra, ct Kom 
YNOs-N 
2 4 
mg/l 0.0] 2.0] 4.0| 6.0 


8.0] 10.0] 12.0] 14.0] 16.0| 18.0] 20.0] 22.0] 24.0] 26.0] 28.0] 30.0] 32.0 2 


1.9] 3.9] 5-9] 7-9] 9.9] 11.9] 13.9) 15.9] 17.9] 19.9] 21.9] 23.9] 25.9] 27-9] 29,9) 31-9] 33.9 


0.00—0.09| 4 8 2 3 6 2 | — 2 | — | — I I I 30 
0.10—0.19| 10 9 7 8 5 2 3 2 I 2 | — I | — | — 1 | — | — SI 
0.20—0.29] 5 5 I 2 2 15 
0.30—0.39| 3 5 | — I 10 
0.40—0.49| I I 2 | — I 5 
0.50—0.59| I I I 3 
0.60—0.69| 2 2 
0.70—I.09 
T.TIO—1.19 I I 
2; 27010297 142391274. 7.15 herp acl 206 2 —|—|: —|: A 


total nitrogen amount supplied in kg per ha (s) 
and the amount of precipitation P up to the 
moment to which s refers, the simplest ex- 
pression satisfying the general conditions indi- 
cated by the frequency tables seems to be: 


x= a logn (EP) (1) 


where « is a certain constant, the physical 
meaning of which will be clear from the 
following. 

For the mean concentration in milligrams 
per liter of the substance (NH,-N or 


NO,-N) in the precipitation P, we get 


1 P 
sh neu) (2) 

The concentration o in the falling precipita- 
tion at the moment, when the amount N has 


fallen is given by: 


ds 400" a 
dD. AP >) 


05 21005 


o is here the concentration of the special 
ions in the precipitation falling at any given 
moment. As is seen from (3) o decreases rather 
rapidly with an increase of P. Thus, for in- 
stance when ı mm rain has fallen, o has only 
half of the value which it had just at the 
beginning of the rain. 

The observations available to us are not 
sufficient to justify more detailed conclusions 
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regarding the mathematical expression used. 
It seems probable, however, that the rate of 
change of o with P must to some extent 
depend upon which kind of ions we consider, 
viz. that the effectivity of the rain in “‘washing 
out” various substances from the atmosphere 
may be dependent upon the nature of the 
substance in question. This may be taken 
into consideration by adding a constant and 


by giving (3) the form: 


100 & 
un (4) 
where k may then be a constant characteristic 
of the type of ions in question. 

It would probably be still more satisfactory 
to give o an exponential expression: 


ds 


AP = xe BP (s) 


09 = 
which for s would give: 


as Fal Lee pea DC) 


Here « is as before the concentration at P = 0, 
and ß is a measure of the rate at which the 
concentration decreases with an increase in P. 

Within the interval of P available through 
the observations, either set of formulas, (1) 
and (3) on the one hand and (5) and (6) on 
the other, may be applied. The accuracy of the 
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Table 3. Total amount of NH,-N supplied to the ground (kg/10*m°) as a function of total 
amount of precipitation. All air masses 


SNHi-N 
kg/ha 
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Total amount of NO,-N supplied to the ground in kg/Io!m? as a function of total 
amount of precipitation. All air masses 


0.0 


2.0) 4.0] 6.0] 8.0] 10.0| 12.0] 14.0] 16.0| 18.0| 20.0 


Pire Cp it a tro Rene 


22.0| 24.0| 26.0] 28.0| 30.0] 32.0 


11.9| 13.9] 15.9] 17.9] 19.9] 21.9 


23.9] 25.9] 27.9] 29.9} 31.9} 33.9 


H 


27 


29 


13 


14 


Zu. 


I 


3 


ON THE CONTENT OF NITROGEN IN ATMOSPHERIC PRECIPITATION 35 


results is practically the same. If we pass to the 
upper limit (P = co, or P very large) there 
is, however, a fundamental difference between 
the two equations (1) and (6). Both are alter- 
native expressions for the amount of NH,-N 
or NOs-N conveyed to the soil through 
the rain. In both cases the actual concentration 
: ds 
at a given moment, IP’ 
P increases above all limits. But if using 
formula (1) the quantity of substance con- 
veyed to the soil increases without upper 
limit when P gets larger and larger, while 


approaches zero when 


the expression (6) shows an upper limit, F 
which will never be surpassed. The difference 
has at present more of an academic than a 
practical sense; probably neither of the for- 
mulas discussed will be able to give a true 
expression outside the limited interval cov- 
ered by our observations. Therefore we will 
for the reasons indicated above, base the fol- 
lowing discussion on the simple formulas (1) 
to (3), which seem to reflect rather well the 
main results of the observations on the content 
of NH,-N and NO,-N. 

Applied to the 117 measurements of the 
content of NH,-N, the equation (1) takes 
the form: 


SNEL-N = 0.0275 log (1 + P) (7) 
from which we get: 


log (1 + P) 


EN 8 
YNH«-N = 2-75 P (8) 


For the content of NO,-N, we obtain in a 
similar way: 
SHON ='0.0132ulog—(1 +P) (9) 


log (1 + P) 
= (10) 


YNO-N 1.32 


In Table 5 computed values of concentra- 
tion are given for values of P up to 40 mm. 
From the table we find for instance that the 

"mean concentration (y) corresponding to I 
mm precipitation is about 8 times larger than 
the mean concentration in a rain or snow fall 
of 40 mm. 


Table 5. Computed mean values of nitrogen 
concentration for all air masses as a function of 
total amount of precipitation 


NH,-N NO,-N 
p ee ae 
mm A Ss We s 
mg/l kg/ha mg/l kg/ha 
I 0.83 0.008 0.40 0.004 
3 0.55 0.016 0.26 0.008 
5 0.43 0.021 0.20 0.010 
10 0.2 0.029 0.14 0.014 
15 0.22 0.033 O.II 0.016 
20 0.18 0.036 0.09 0.017 
25 0.16 0.039 0.07 0.019 
30 0.14 0.041 0.06 0.020 
40 O.II 0.044 0.05 0.021 


The total supply to the ground of the ions in 
question covers the period I. 3. 1948—1. 3. 
1949 and 1.8. 1949-31. 12. 1949.: Not all 
precipitation cases, however, were investigat- 
ed in this study. The actually measured quan- 
tities supplied to the ground amount to: 


For NH,-N : 2.37 kg per ha 
» NON E09 2 I» 


The number of cases of precipitation during 
the period mentioned is, however, somewhat 
larger than the number on which measure- 
ments were made, namely 190 against 117. 
In order to find the real supply corresponding 
to all cases, we may use the relation between 
supply and amount of precipitation computed 
above and apply it to the frequency distribution 
of the precipitation within various intervals. 
This frequency distribution is as follows: 


P mm Number of cases 
O.I—T 53 

I—5 82 
Se 36 
10—20 14 
20—40 5 


Total = 190 


Taking the mean values of P within every 
interval as independent variable and using the 
formulas (7) and (9), we obtain through mul- 
tiplication with the number of cases within 
each interval and through adding all cases, 
the following amounts: 
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For NH,-N : 3.22 kg per ha 
» NO,-N :1.549 » >» 


We may now proceed a step further and 
estimate the average supply during a normal 
year. For the present purpose a computation 
based on a period of 10 years seems sufficient. 
Investigating the frequency distribution of 
various cases of precipitation at Uppsala, 
which is situated close to Ultuna, we obtain 
the following table: 


P mm Number of cases 
0.1—I 2166 
I—5§ 2207 
$—10 664 
10—20 302 
20—40 54 
40—60 8 


Using these data in a similar way as above we 
obtain for the average annual supply: 


For NH,-N : 1.90 kg per ha 
» “NO Z-N*: 0.9197 35 


The number of observations treated above is 
not so large as to allow us to attach more 
than a moderate weight to these numbers. 
Several circumstances may work in the direc- 
tion to make the generalization, which we 
have attempted in order to compute average 
values, rather uncertain. It seems rather pro- 
bable that the constant of the formulas (r) 
may be subjected to rather large variations 
from time to time, not only of a seasonal but 
also of a secular character. The fact, which 
will be discussed below, that we undoubtedly 
have a variation of « from one type of air 
mass to another, naturally indicates that long 
periodic changes in the frequency of certain 
air masses will reflect themselves also in 
the values of y, which in their turn are meas- 
ures of the concentration of the ions in the 
precipitation. 


2. Content of NH,-N and NO,-N in precipita- 
tion within different air masses. 


A preliminary investigation of the depen- 
dence of the nitrogen content on the time of 
the year and on the occurrence of thunder- 


storms gave on the whole a negative result. 
This may, however, be due partly to the limi- 
tation of the observations in time and space. 

We then considered it to be of great interest 
to look for a possible connection with the 
type of air mass. A direct investigation of the 
physical conditions of the air masses in which 
precipitation occurs, viz. their temperature 
distribution, the height and character of 
clouds in the case of precipitation would pro- 
bably be still more illuminating as regards the 
character of the physical process which fur- 
nishes the nitrogen content here in question. 
However, in the absence of accurate and 
sufficiently detailed data concerning these 
factors at the place of observation, we shall at 
present confine ourselves to determine whether 
or not the influence of various physical factors 
may be reflected in connection with the rough 
synthesis of their influences, expressed through 
the air mass type. 

In designating the type of air mass in the 
117 cases of precipitation we had the very 
valuable help of meteorologists of the fore- 
casting section of the Swedish Meteorological 
and Hydrological Institute. Especially to Mr 
Bengt Magnusson, chief forecaster, we have 
the pleasure of expressing our most sincere 
thanks. 

A very detailed indication of air mass is 
scarcely justified by the material as it is. A divi- 
sion into three chief types, namely tropical, 
polar and arctic air, seemed to meet the require- 
ments of a preliminary analysis. The result 
was that in the 117 cases of precipitation, 20 
took place in tropical, 92 in polar and 5 in 
arctic air. Whether this proportion ought to 
be considered as normal is difficult to say, but 
there is no reason to assume that it deviates 
considerably from average conditions. In 
Tables 6—9 the frequency distribution of the 
concentration is shown with its dependence 
on air mass in the case of the two nitro- 
gen ions NH} and NO, The arctic air 
which is represented by only 5 cases is not 
included in the tables. Tables 10—13 show the 
corresponding frequency distribution for the 
supply (s) of nitrogen. 

We find in all these cases the same diagonalic 
distribution as was previously discovered when 
the precipitation was studied independently 
of air mass. 

A computation of the coefficient « gives for 
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Table 6. Frequency distribution of the mean concentration of NH,-N (mg/l) in tropical air masses 
given as a function of the total amount of precipitation 
a Ee oe a De nn © 
PRE Cp ta tion ment 
YNH«-N 


mg/l 0:02 122,02 124.0317770:0 1.8.0, | 10:01 12.0]f 74:01 16.0 \|/ 18.0 | 20:0 || 22.0 2 


291331259) 7.9 | 9.9 | 12.9] 13.9) 15.9) 17.9.| 10.9 | 21.9 |. 23.9 


0.00—0.19 I = == 2 B 
0.20—0.39 — I I == — I I I 5 
0.40—0.59 I — == = I 2 
0.60—0.79 I I a — is 3 
0.80—0.99 — 2 2 
1.00—1.19 — 3 B 
1.20—I.39 I — I 2 

2 | 3 | 7 | 2 — | 3 | I I 2 | | | I | 20 


Table 7. Frequency distribution of the mean concentration of NH,-N (mg/l) in polar air masses 
given as a function of the total amount of precipitation 


Bre Chpt ation 


YNHi-N : a ee x 
en 0.0] 2.0] 4.0] 6.0] 8.0] 10.0] 12.0] 14.0] 16.0] 18.0] 20.0] 22.0] 24.0] 26.0] 28.0] 30.0] 32.0 
D 
1.9] 3-9] 5-9} 7-9] 9.9/11.9| 13.9] 15.9] 17.9] 19.9) 21.9] 23.9] 25.9] 27.9] 29.9] 31.9) 33.9 
0.00I— 
2 2 | — — 2 — | — — 28 
en 4 3 2 6 4 2 I I I I I 2 
0.20I— 
8 5 5 = 2 2, 
0.391 zu SA a ae 23 3! I 3 
0.40I— AIS ss ns 
0.591 3 5 2 
0.601— > 2 5 hi 
0.791 
0.80I— 5 3 AE i 6 
0.991 
1.00I— 
I I 
1.191 
1.20I— à 
1.391 
1407 |) 1 : 
1.591 
1.001 rs. 4 5 
1.791 
2} 224 6270 gene 0232 62 2 2 2 I 2 I I —|—-|: — | : 92 


the different air masses in the case of NH,-N For NO,-N we get: (cf. figs. 3 and 4): 
(cf. figs. ı and 2): 


Tropical air : « = 0.0354 Tropical air :« = 0.0142 
Polar PRIE 10 0202 Polar D 0,0131 
Arctic ee == O10172 Arctic » 1% = 0.0094 
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Table 8. Frequency distribution of the mean concentration of NO,-N (mg/l) in tropical air masses 
given as a function of the total amount of precipitation 
a ee ee EE SS SS er — 


Pirecipita tom mm 
} nt 
NEN 0:0 | 2.0 | 4.0 | 6.0 | 8.0 | 10.0] 12.0] 14.0] 16.0 | 18.0 | 20.0 | 22.0 Z 
mg/l : 
1294 03:9, || Gaon | Ke) |] KOH) nr 0 trs 0 ES 0 77:94 5279-97 TON RS 
0.00—0.09 — I — | — — | — 2 4 
0.10—0.19 — I I — I I I 6 
0.20—0.29 I —| —| — 2 
0.30—0.39 I 4 5 
0.40—0.49 — I I 2 
0.50—0.59 = 
0.60—0.69 I a 
SRE 
Table 9. Frequency distribution of the mean concentration of NO.-N (mg/l) in polar air masses 


given as a function of the total amount of precipitation 


Patera py i tra. temoin 


YNOs-N 
mg/l 


mm 


0.0] 2.0] 4.0] 6.0} 8.0| 10.0] 12.0| 14.0| 16.0] 18.0] 20.0] 22.0] 24.0] 26.0] 28.0] 30.0} 32.0 x 


I.9} 3.9 


5-9] 7-9} 9-9) 11.9 13175-977979 EAN FEIN Sn 15.9| 17.9| 19.9] 21.9| 23.9] 25.9| 27.9| 29.9] 31.9| 33-9 


aN 


tv 
H 
+ 
On 


Ze | elle LS 2 2 2 2 


As « can be regarded as a measure of the 
average content reduced to a certain quantity 
P of precipitation, it is evident that the 
tropical air in general contains more nitrogen 
ions than polar air, and considerably more 
than arctic air. If we denote « for the tr opical 
air with 1.00, the values for polar and arctic 
air are the following: 


Air mass NH,-N NO,- 
Tropicals. 7 1.00 1.00 
POlAt enr 0.74 0.92 
ATCUCRE RE 0.49 0.66 


3 6 2 3 5 I I I 23 
10 8 6 8 I | 
4 5 I I I 12 
2 I | — I I _ 5 
PSN MENT 3 
I I I = : é 3 
ee ee EN es : 


| 


I 


; |: 


— | : 


The table is of interest in throwing some 
light upon the i importance of various air masses 
for the supply of nitrogen. These conditions 
are of great interest soie it comes to evaluating 
the supply of nitrogen at places with different 
geographical position; for instance for making 

clear the differences between the northern and 
southern parts of Sweden. On the other hand 
the difference seems in that case to be larger 
than would be implied from our present in- 
vestigation. Northern Sweden receives less 
tropical air than southern Sweden, which 
also is true in case of precipitation. The supply 
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Table 10. Total amount of NH,-N supplied to the ground (kg/ro!m?) as a function of total 
amount of precipitation. Tropical air masses 
Breet pitation mm 

SNHi-N 

kg/ha 0:05 02:0 1674500) 36:0 5188.0 |7.10:0| 212,01 1440) 16.0 || 18.0 20101127 0 x 
7,30 02:92 RE Os 16295 NO OL TON ELS OF 15.9 1777.79 19.071" 21g 23-9 

0.00I—0.010 I I 2 
O.01I—0.020 2 I I — I 5 
0.021—0.030 | — 3 I 4 
0.03I—0.040 = 2 u i I — I 5 
0.041—0.050 - I I 
0.05 1—0.060 —- — I — I I 3 
eee ame we le 
Table 11. Total amount of NH,-N supplied to the ground (kg/10‘m?) as a function of total 


amount of precipitation. Polar air masses. 


Binecipıtatıon 


6.0] 8.0} 10.0] 12.0] 14.0 


mm 


16.0] 18.0] 20.0] 22.0] 24.0] 26.0] 28.0] 30.0] 32.0 a 


7-9] 9.9) 11.9| 13.9) 15.9 


17.9] 19.9] 21.9| 23.9] 25.9] 27-9] 29.9] 31.9] 33.9 


30 


15 


Se ale 


of nitrogen through precipitation is on that 
account larger in the south than in the north. 
But the difference is more pronounced, for 
the tropical air arriving in the northern parts 
seems to have a lower nitrogen content than it 
has at a lower latitude; whether from the 
effect of being diluted with polar air or from 
being already somewhat washed out by preci- 
pitation during its journey is difficult to say. 
Probably both these causes are at work. 
Finally, if we compute the ratio between 


the deviation of the observed values from the 


computed curve, according to formulas 1—3 
on the one hand, and the standard deviation 
of the observed values on the other, we get 
the values of the following table: 


e 


Air mass NH,-N NO;-N 
Tropwalsee. set 2% 1.27 
Polar sean uit 2.028 re 
ATEN este 1.03 0.95 


With exception for the arctic air with its 
weak representation in the material, the rela- 
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Table 12. 


Total amount of NO,-N supplied to the ground (kg/rom?) as a function of total 


amount of precipitation. Tropical air masses 
ed 


Precıpitatrom 


mm 


oa eS SSS EEE 


SNO:-N 
kg/ha 


2 


14.0| 16.0 | 18.0 | 20.0 | 22.0 


11.9| 13.9 


1591 £729 | 19.9 |) 27.907 23.9 


0.00I—0.005 2 


0.006—0.0IO I 
0.01I—0.015 u 


0.016— 0.020 — 
0.02I—0.025 


HHWN 
nln | 


| 
H 
HWE OWN 


0.026—0.030 


0.031—0.035 


Bere 


Table 13. 


| I 20 


Total amount of NO,-N supplied to the ground (kg/10m°) as a function of the total 


amount of precipitation. Polar air masses 


Precip ita ti Om) mm: 


0.0] 2.0] 4.0] 6.0] 8.0] 10.0] 12.0] 14.0] 16.0] 18.0] 20.0] 22.0] 24.0] 26.0 28.0 


| 30-0] 32.0 x 


1.9] 3-9] 5-9} 7-9] 9.9} I1.9| 13.9) 15.9] 17.9] 19.9] 21.9] 23.9] 25.9] 27.9 29.9} 31.9 33-9 
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tion of these numbers to one another seems 
to form a support for the relevancy of the 
mathematical treatment. 

The figures included here graphically give a 
comprehensive idea of the content of the 
tables and the mathematical expressions used. 


II. Summary of results and survey of some 
unsolved problems 


In the present paper we have shown: 
(1) The amounts of NH,-N and NO,-N 


contained in precipitation are dependent upon 
the amount of precipitation which has pre- 
viously fallen. 


As the precipitation continues to fall it will 
contain fewer and fewer of these ions. 


(2) A mathematical expression is given 
for the dependence of concentration on the 
amount of precipitation fallen. 


(3) With aid of this relation a computation 
has been made of the normal supply of 
NH,-N and NO;-N at Uppsala (Ultuna) 


> 


CPS 
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0.07 =. 
NH,—N 
0.06 
0.05 } 
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0.04 
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Fig. 1. Content of NH,-N in precipitation within 
different air masses. 
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A= Arctic air 
04 
7 
Pp 
02 A 
mg/y 
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Fig. 2. Mean concentration of NH,-N in precipitation 


within different air masses. 


using known normal frequencies of different 
amounts of precipitation. 

(4) The normal annual supply at Uppsala 
is thus found to be: for NH,-N : 1.90 kg/ha; 
for NO,-N : 0.91 kg/ha. 

(s) A pronounced dependence of the con- 
centration on the character of the air mass has 
been found to exist. Thus the tropical air 
contains between 10 and 30 per cent greater 
amount of the nitrogen ions of the investigated 
kind than the polar air and almost double the 
amount found in the arctic air. 


0.030 = aS er 


0.025 en: = a 


0.020 


0.015} 


0.010 


= Tropical air 
0.005 = > P= Polan.ain Je 4I 
A= Arctic air 


0 mm 5 io 15 20 275 CONS 


Fig. 3. Content of NO;-N in precipitation within 


different air masses. 


ITR 


T - Tropical air 


P - Polar air 


| A- Arc 


mg/v| x =I 
0 mms 10 15 20 25 30 35 


Fig. 4. Mean concentration of NO;-N in precipitation 
within different air masses. 


The investigation reported on here is mainly 
of statistical character. No attempt has been 
made to go into the ultimate causes of the 
relations which have been established. Are 
the nitrogen ions treated here bound to some 
kind of solid or liquid particles on which the 
water particles perhaps condense (nuclei) or 
are they present to a considerable extent in a 
gascous state in the atmosphere and absorbed 
by the water particles to which the affinity cer- 
tainly must be great? No definite answer to 
this question can be given here. On the other 
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Table 14. Computed mean values of nitrogen 
concentration for different air masses as a function 
of total amount of precipitation 


Air P en 
es | agi msi | gba | men kg/ha | mg/l ES kg/ha 
I 10/2100 770.43 0.004 
4 3 0.70 |7,0:02%. 50.28 0.008 
2 5 0.55 | 0.028 | 0.22 0.011 
= 10 0371100 037 120.75 0.015 
Yar 15 0.28 | 0.043 | O.I1 0.017 
= 20 0.23 | 0.047 | 0.09 0.019 
= 25 0.20 | 0.050 | 0,08 0.020 
30 0.18 | 0.053 | 0.07 0.021 
ee ET N 0.059 | 0.06 0.023 
I 0.79 | 0.008 | 0.39 0.004 
3 0.53 | 0.016 | 0.26 0.008 
4 5 0.41 | 0.020 0.20 | 0.010 
‘8 10 0.27 | 0.027 0.14 | 0.014 
tal 15 0.21 0.032 0.11 0.016 
© 20 0.17 | 0.035 0.09 ! 0.017 
A 25 0.15 | 0.037 0.07 | 0.019 
30 0.13 | 0.039 0.06 | 0.020 
40 0.11 ee 0.05 | 0.021 
I 0.52 | 0.005 0.28 | 0.003 
3 0.35 0.010 0.19 0.006 
_ 5 0.27 10. 0.013 0.15 0.007 

Re) 

& Io 0.18 | 0.018 0.10 | 0.010 
aS 15 0.14 | 0.021 0.08 | 0.011 
5 20 0.11 | 0,023 0.06 | 0.012 
a 25 0.10 | 0.024 0.05 | 0.013 
30 0.09 | 0.026 0.05 0.014 
40 0.07 | 0.028 0.04 | 0.015 


hand the almost constant ratio of about 2:1, in 
which the NH,-N, and NO,-N respectively 
appear in the precipitation seems to suggest a 


common origin — perhaps some kind of nuclei, 
in which they enter in an almost fixed ratio. 


How is the concentration at the end of 
very heavy rain or snow falls? Does the con- 
centration approach zero or a negligeable 
value, or will there arise some kind of equi- 
librium conditions, wherein as many nitrogen 
ions of the investigated types are produce 
through some process, as disappear through 
the “washing out” of the precipitation? 


What is the reason for the greater concen- 
tration in tropical air masses than in polar 
and arctic ones? The facts seem to suggest a 
supply from the earths surface in the same 
way as pollen and dust and spores are brought 
to the air through the stronger convection 
processes at low latitudes. It seems also pos- 
sible that the existence of nitrogen ions in 
the atmosphere is to some extent connected 
with electrical phenomena, as thunderstorms 
and other forms of discharges, which occur 
more frequently at low than at high latitudes. 
No definite conclusions can, however, be 
reached from the limited material treated 
here. 


An investigation which suggests itself is on 
the connection between certain biological and 
agricultural phenomena, which are dependent 
upon the nitrogen supply as for instance plant 
growth and ion and on the other 
hand the occurrence of precipitation from var- 
ious kinds of air masses as defined through 
their frequency distribution. Such an investiga- 
tion falls, however, outside the scope of our 
present field of ach. 
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The Distribution of Temperature and Wind Connected 

with Active Tropical Air in the Higher Troposphere 

and Some Remarks Concerning Clear. Air Turbulence at 
High Altitude 
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Abstract 


A study is made of the distribution of wind and temperature in the higher troposphere and 
the lower stratosphere in cases with active Tropical air in middle latitudes. Only observed 
winds are used. A model for extending the fronts into the lower stratosphere is given. 

In connection with the proposed model a case with observed clear air turbulence at high 
altitude is studied. Some remarks on the statistical treatment of the reported cases with clear 


air turbulence is made. 


The concept of fronts has been used in syn- 
optic meteorology for about 30 years as far as 
sea level charts are concerned. In synoptic 
aerology, however, the fronts have not been 
fully accepted, although J. Bjerknes, Palmén, 
Van Mieghem, and others already in the 30's 
proved that is was possible to use this concept 
also in the middle troposphere. The aerological 
work (in Europe) was founded, however, to a 
great extent on aircraft ascents, the ceiling of 
which was about $ km. It was as a rule not 
difficult to identify the fronts up to this height, 
but when the soundings began to reach 
heights between 10 and 20 km new problems 
arose. It was firstly the swarm ascents, made 
possible by the Jaumotte meteorograph, and 
secondly the appearance of the radiosond that 
brought to the fore the problem of prolonga- 
tion of the fronts into the higher troposphere 
(and the lower stratosphere). 


1 Written while on leave of absence from Swedish 
Meteorological and Hydrological Institute, Stockholm. 


It is especially when studying the three- 
dimensional structure of the atmosphere by 
means of vertical cross sections one meets this 
special problem. 

The rapidly increasing importance of syn- 
optic aerology during World War II, brought 
about especially by the development of air- 
craft efficiency, focussed the meteorologists’ 
attention still more on conditions in the higher 
troposphere and the lower stratosphere. 

This paper is a preliminary report of an 
investigation concerning the problem of ex- 
tending the fronts into the lower stratosphere. 

This problem has been treated by J. BJErK- 
NES and PALMÉN (1937), J. BJERKNES, Mirp- 
NER, PALMÉN, and WEICKMANN (1939), Ny- 
BERG and PALMÉN (1942), but the available 
material was mostly sparse, and wind observa- 
tions from the critical region were entirely 
missing. The present author has taken up the 
problem because of the fact that the radiosond 
network nowadays by far exceeds that, on 
which the earlier works were founded. 
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„ TROPICAL AIR 


POLAR AIR 


Fig. 1. Two soundings representative for Tropical and 
Polar air, respectively, in western Europe and over the 
middle Atlantic in November. 


One week of November 1949 is treated by 
the author for another purpose. From this 
material two soundings were constructed?, one 
representing Tropical air and the other Polar 
air (see fig. 1). The temperature difference 
between the two air masses’ increases with 
height to about 400 mb, from 7° C to about 
16° C. Above the tropopause (at 325 mb) the 
Polar air is isothermal, and at 257 mb both air 
masses have the temperature of — 49.0° C. 
Then the horizontal temperature gradient is 


2 The details of the computations will be accounted 
for in the complete report, which will be published in 
Arkiv for Geofysik, vol. 2 (in print). 


reversed and increases up to 130 mb, where 
the difference is about 11° C. Now, the question 
is whether this stratospheric temperature dif- 
ference is concentrated to a front, as is the case 
in the troposphere, or not. Furthermore it 1s 
interesting to study the behaviour of the front 
in the upper part of the troposphere. 

In order to study these phenomena I have 
chosen two cases, one from November 1949, 
and another from October 1949. 

On the 7th of November 1949, an almost 
stationary front was located over western 
Europe, extending towards westnorthwest 
to the middle of the Atlantic. This front began 
to move northeastwards over the Atlantic, 
later towards eastnortheast. A rather strong 
sea-level low developed during this period, 
and on the oth at 0300 GCT its warm front 
had just reached northern Ireland and western 
France (see fig. 2 a). The cold front extended 
from westernmost Ireland down to the Azores. 

A high over Russia was almost stationary, 
and the southeasterly to southerly air current 
at its southwestern and western side was very 
stable and had persisted for about four days. 

The cold air over Central Europe between 
these two warm currents was maritime Polar 
air, which had come from the Atlantic and 
which had partly been Arctic air. 

On the soo mb chart (fig. 2b) the advection 
of the Tropical air over and east of the British 
Isles is very pronounced. The warm front is 
especially well-defined over northwestern 
Germany, the Netherlands, and Belgium, 
and here the density of the radiosond network 
is so great that the analysis may be considered 
as quite reliable. 

The broad southerly current over eastern 
Europe is well developed, but the front at 
500 mb connected with this current is not as 
strong as the western one. 

As a whole the sea level chart as well as the 
500 mb chart in fig. 2 gives a good picture of 
an upper-air pattern, which is rather common, 
at least if western and eastern Europe are con- 
sidered separately. 

In order to illustrate the distribution of wind 
and temperature in this case, I have drawn a 
vertical cross section parallel to the latitude 
circle of about 54.5° N. The section extends 
from Valentia in western Ireland to the British 
Zone in Germany. The cross section can not 
be extended further eastwards, because the 
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Fig. 2a. Sea-level chart for 0300 GCT 9 November 1940. 
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Fig. 2 b. soo mb chart for 0300 GCT 9 November 1949. Solid lines are isohypses for every 4 geodynamic 
decameters, dashed lines are isotherms for every 2° C. 
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Fig. 3 a. Ascent curves, 0200 or 0300 GCT 9 November 
1949, for Hannover and Iserlohn (British Zone in Ger- 
many), Uccle (Belgium), De Bilt (the Netherlands), 
Downham Market, Larkhill, and Valentia (British Isles). 
The left curves are the ascent curves for Hannover 
and Jever (British Zone in Germany) for 2000 GCT 
8 November, and the second from the right is the 0900 
GCT ascent from Larkhill on the 9th of November. 


9 x! 1949 
03° GCT 


Fig. 3 b. Vertical cross section from Valentia (Ireland) to 
Hannover (British Zone in Germany) for 0300 GCT 
9 November 1949, approximately along 54.5°N. Thin 
dashed lines are isotherms (° C), thin solid lines isovels 
for observed wind speed (m sec~!, projected on a 
direction of 330°). Intermediate values of temperature 
and wind speed are given by dotted and dashdotted lines, 
respectively. Heavy lines represent tropopauses, frontal 
boundaries, and inversions. Thin vertical lines give the 
location and vertical extent of temperature soundings.’ 


Fig. 3 c (next page). Vertical cross section as in fig. 3 b, 
giving isentropes (thin dashed lines, °K) and isovels 
for observed wind speed (thin solid lines, m sec-1). Inter- 
mediate values for the wind speed are given by dash- 
dotted lines. Frontal boundaries, tropopauses, and inver- 
sions given as heavy lines. Thin vertical lines give the 
location and vertical extent of the wind soundings. 
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frontal system of eastern Europe would then 
enter the picture. To avoid this difficulty we 
have to some extent changed our space cross 
section to a time cross section in the eastern- 
most parts. The changes within the Polar air 
during six hours are quite small and it is there- 
fore permissible to use the Hannover and Jever 
soundings from the 8th of November, 2000 
GCT, in our section. The ascents are plotted 
so as to be in the right positions vis-à-vis the 
front, the velocity of which parallel to the cross 
section was observed to be about 35 km ht. 

Furthermore the sounding from Penzance 
at 0300 GCT on the oth of November is 
a little too warm in the higher troposphere. 
This fact has been revealed through the analy- 
sis of the upper-air charts, and we have there- 
fore substituted the Penzance ascent by the one 
from Larkhill, at 0900 GCT on the same day. 
If we use the estimated velocity of the front 


the latter ascent will be situated approximately 
at the right place. 

The ascent curves of this cross section are 
given in fig. 3 a. The Polar front is easily seen 
in all soundings except those three to the left. 
The inversion in the Hannover ascent from 
0300 GCT may be taken as an indication of 
the front, but the earlier history of this air as 
well as the wind distribution given by this 
ascent, does not support this solution. This 
question will be subject for a more detailed 
discussion later in this paper. 

The tropospheric front is very distinct in the 
isotherm field (see fig. 3 b). It 1s a particularly 
interesting fact that practically all soundings 
give inversion in the lower portion of the frontal 
layer and weak vertical temperature gradient 
in the upper portion. In the fully developed 
stage almost all soundings through this warm 
front give the same picture. Studying the ascents 
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from the 8th through the same warm front 
one comes to the conclusion that this fact may 
be taken as an indication that the vertical 
motions on one side of the front are rather 
strong. A discussion of this interesting question 
would, however, lead too far. 

The tropospheric Tropical air is nearly baro- 
tropic, except in the neighbourhood of the 
front and in the westernmost part, where the 
prefrontal cooling ahead of the cold front can 
be seen. The Polar air is slightly barocline in 
the troposphere and near the upper parts of the 
tropospheric front it is quite strongly barocline. 

In the stratosphere the isobaric temperature 
differences within the Polar and the Tropical 
air, respectively, are rather small, except in the 
Tropical air near the front. The contrast be- 
tween the two air masses amounts to about 5° C 
and seems to be concentrated to the front. 

Fig. 1 showed the well-known fact that the 
Tropical air is colder than the Polar air above, 
say, 350—250 mb. If the Tropical air is active 
throughout the troposphere and lower stratos- 
phere, then the Polar front must act as a warm 
front in the troposphere and as a cold front 
above the lower tropopause in the Tropical air. 

Throughout this paper we have used the 
expression ‘stratospheric Tropical air’ and the 
like. According to our opinion it is not correct 
to regard the tropopause as the upper limit of 
the Tropical and Polar air (in middle and high 
latitude). The tropopause (or tropopauses) is (are) 
mostly a phenomenon (phenomena) occurring within 
the air mass in question. 

The wind field (figs. 3 b and 3 c) contains a 
strong jet at about the 350 mb level with a wind 
speed of about 75 m sec~!. The figures give 
the component of the wind speed parallel to 
330°, i.e. a direction almost perpendicular to 
the cross section. The vertical wind shear 
through the tropospheric front is about 13 m 
sec"! km“! in the lower part and about 33 m 
sec-!km~! between the 400 and soo mb levels. 

The isobaric wind shear is considerable. The 
anticyclonic shear is about 9 m sec-1 (100 
km)-t, if we assume a linear fall of the wind 
speed from the jet axis to the station Downham 
Market (DM). The real shear is, however, 
probably greater in a region somewhere in 
between, where values approaching the Corio- 
lis parameter are likely. 

The cyclonic wind shear is definitely greater. 
Assuming a linear velocity field between the 
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stations Iserlohn (IS) and Hannover (HA) we 
get 32 m sec! (100 km)-!. Comparing the 
velocity distribution both in horizontal and 
vertical direction. where the Iserlohn ascent 
passes through the warm front, we find justi- 
fication for concentrating most of the isovels 
to the front. In this manner we obtain an 
isobaric cyclonic wind shear of about 60 m 
sec-! (100 km)-}, or.more than five times the 
Coriolis parameter. 

There is, unfortunately, one complication 
that ought to be discussed a little. The launch- 
ing time of the ascents may deviate somewhat 
from what is indicated in the synoptic reports, 
and the observed wind speed has both system- 
atic and accidental errors. Assuming a time 
error of 30 minutes and a speed error of Io m 
sec"! and using the estimated velocity of the 
front, we obtain 22 m sec-1 (100 km)-1 as 
value for the cyclonic wind shear, with a linear 
fall of the wind speed. The assumed errors are 
certainly too great in most cases and therefore 
we may give the interval 45—60 m sec-! 
(100 km)! for the wind shear as a convenient 
one in this case. 

The isobaric wind shear is thus greatest in 
the region where the temperature difference 
disappears and the front becomes vertical. Thus, 
with the aid of the wind observations it should 
be possible to trace the front between 400 and 
250 mb. In this region it is often difficult to 
analyse with the aid of temperature alone, 
because the temperature distribution is blurred 
by strong vertical motions. (In the complete 
report we are going to discuss to some extent 
the vertical motions connected with this 
model.) 

Here it must be emphasized that throughout 
this paper we use the ordinary definition of a 
front as a discontinuity in density. It is only 
in the region where the temperature difference 
between the two air masses disappears that 
we use a special definition of the front. In this 
certainly rather thin layer we define the front 
as the region with great horizontal cyclonic 
shear. 

As item two we have chosen a case from 
October 1949. During the period 22nd to 24th 
unstable, fresh maritime Polar air covered the 
greater part of northwestern Europe. Between 
Greenland and the Scandinavian peninsula 
there was a continuous northerly flow of Arctic 
air, which was transformed on its way south- 
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wards over the ocean. It was therefore not 


possible to trace an Arctic front in this region. 
Over the Atlantic Ocean south of so °N a 
series of disturbances moved eastwards to the 


southern parts of the British Isles and France. 


One of these systems which on the morning 
of October 25th was situated south and south- 
west of the British Isles, became more intense 
than the preceding ones. The warm front of 
this system was very distinct and can easily be 
found on the sea-level chart. (The sea-level 
fronts are indicated in fig: 4 with the interna- 
tional symbols). 

The 500 mb chart (fig. 4) shows a strong 
southwesterly current over southwestern Eng- 
land and adjacent seas. This current veers west 
and passes over northern Germany, for instance, 
as a due westerly wind. The isotherms are 
crowded where the Polar front intersects the 
5oo mb level, and the temperature difference 
between the Tropical and the Polar air through 
the front is about 10° C. 

This time I prefer to show a time cross sec- 
tion instead of a space cross section. The Brit- 
ish radiosond network, with its six hours 
intervals, is well fitted for this purpose. In fig. 
5 b I have given a time cross section for Alder- 
grove, the radiosond station in northern Ire- 
land, and for the period 0300 GCT 24th—0900 
GCT 25th of October 1949. 

The ascent curves given in fig. 5a have 
essentially the same characteristic features as 
have those of fig. 3 a. The warm front is easily 
detectable from 1500 GCT on the 24th and 
thereafter. The deep layer between 418 and 
330 mb in the 0900 GCT ascent is not be re- 
garded as the frontal layer. The lowest charac- 
teristic point (418 mb) must be identified with 
the lower tropopause of the soundings six 
hours earlier. 

The isotherm field of fig. sb gives as a 
whole the same picture as the one in fig. 3 b. 
As in the first case, the difference in vertical 
stability on both sides of the front between, 
say, 350 mb and 200 mb is very striking. Com- 
paring the ascents from o900 GCT and 2100 
GCT on the 24th, we get the following figures 
for the vertical gradient of temperature an 
potential temperature. The ascents give (the 
figures for Tropical air are given within 
brackets) a temperature difference of — 8.9° C 
(— 23.9°C) and a potential temperature dif- 
ference of + 41.0° C (+ 18.5° C). The mean 
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Fig. 4. soo mb chart for 0300 GCT 25 October 1949. 

Isohypses are given as solid lines for every 4 geopotential 

decameters, isotherms as dashed lines and for every 2° C. 

The fronts at the sea-level are given with international 
symbols. 


vertical distance between the 200 and the 350 
mb levels is 3,713 m. Thus, we get as values 
for the vertical gradient of temperature and 
potential temperature — 2.4° C km! (—6.4°C 
km) and, ++ 11.0° Ckm=! (+ 5.0° C km=t), 
respectively. ; 

The difference in wind velocity on both 
sides of the front between 300 and 250 mb 
is great. In the Polar air just ahead of the front, 
there is a wind speed of about 40 m sec-1. In 
the Tropical air on the other hand the jet seems 
to reach a value of at least 80 m sec. 

If we now take into consideration (1) the 
great differences in vertical stability, (2) the 
observed wind speed on both sides of the front 
and (3) the difference in temperature between 
the two air masses, measured along an isobar, 
below and above the height of zero tempera- 
ture difference (in fig. 3 b at about 325 mb and 
in fig. 5 b at about 310 mb) we reach the fol- 
lowing conclusion. There is no reason not to 
extend the tropospheric warm front into the lower 
stratosphere, where it acts as a cold front. Between 
about 350 and 200 mb the front is a boundary 
between tropospheric Tropical air and stratospheric 
Polar air. 


Remarks on clear air turbulence 


Finally I want to discuss a special pheno- 
menon, which in many cases seems to be con- 
nected in one way or another with the atmos- 
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Fig. 5 a. Ascent curves for the British radiosond station 
Aldergrove, northern Ireland, for the period 0300 GCT 
24 October—1s500 GCT 25 October 1949. 


pheric model discussed above. During the 
last ten years clear air turbulence above 5 km 
has been subject of investigations, especially in 
Great Britain and U.S.A. We shall start with 
a brief summary of two papers, one by Histor 
(1951) and the other by BANNON (1951), both 
summarizing the British investigations. 

A series of experimental flights were made 
by British aeroplanes between 6 and 12 km 
above sea level in order to try to get quantita- 
tive data of the clear air turbulence.! This type 
of turbulence? was not infrequently observed 
over the British Isles above a height of 6 km. 


1 In these investigations only gusts giving a vertical 
acceleration > 0.2 g are considered. The type of air- 
craft used determines to some extent what kind of tur- 
bulence element will be found in the atmosphere. Spit- 
fires and Mosquitos were used and the turbulence eddies 
causing bumps have the dimensions of 15 to 500 meters 
according to Bannon. : 


2 In the British investigations, which were devoted to 
frictional or dynamic turbulence above 6 km, the ther- 
mally induced clear air turbulence in lower layers is 
disregarded. 
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The turbulent layers were usually rather shal- 
low; out of 290 cases 83 °% were observed in 
layers of less than 900 meters depth, but in 8 % 
of the cases the layer was reported to be more 
than 1,800 meters thick. 

The distribution with height shows no great 
variation (58 cases in the layer between 6 and 
7.2 km, 72 between 8.4 and 9.6 km, and ss 
between 10.8 and 12 km). If the material is 
grouped with respect to the height relative 
to the tropopause, then we get a significant 
decrease of frequency above the tropopause 
level. 

The jet streams and their margins are the 
site of many of the reported turbulence cases. 
Bannon tries to combine the Richardson num- 
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Fig. sb. Time cross section for Aldergrove for the 
period 0300 GCT 24 October—0900 GCT 25 October 
1949. Isotherms are given as dashed lines (°C), and ob- 
served winds are plotted in the conventional manner. 
Full barb 5 m sec-1, half barb 2.5 m sec-1, triangular 
barb 25 m sect. Intermediate values of temperature are 
given as dotted lines. Heavy lines are frontal boundaries 
and tropopauses. 
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ber! and the occurrence of turbulence. He 
comes to the conclusion that the bumpiness is 
related to the Richardson number, but the cor- 
relation is not too good. The same condition 
seems to be valid also for the horizontal wind 
shear? contra the frequency of clear air tur- 
bulence. Bannon also concludes that these two 
possible causes of turbulence are to some extent 
independent. 

In an unpublished report by Harrison 
(1950), quoted by Histor, the suggestion is 
made that the wind shear in the horizontal is a 
significant factor. This suggestion is founded 
upon experiences from United Airlines flights 
with DC-6’s between 4.5 and 7.5 km height in 
eR SHAG 

In Bannon’s paper nothing specific is said 
about the connection between clear air turbu- 
lence and the fronts. He reviews a few examples 
of flights during which bumpiness was en- 
countered in inversions or isothermal layers. 

After this short review of the British work 
in this field a synoptic example will be 
given. On April 25th 1949, at 1200 and 
1220 GCT an aeroplane encountered clear air 
turbulence over southeast England. The cross 
section in fig. 6 from Aldergrove in northern 
Ireland to Downham Market in southeast 
England was located quite near the critical 
region. As there are no radiosond measure- 
ments available for 1200 GCT, the 3900 GCT 
soundings are chosen. This cross section gives 
the same picture as do the earlier cross sections 
(cf. figs. 3 b and c, figs. sb and c)?. The front is 
easily traceable and the temperature difference 
between the two air masses is of the order 10°C. 
The cold air in the troposphere is almost baro- 
tropic, except near the frontal boundary, and 
its tropopause is well defined. 
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dz 


1 The Richardson number 


(where g is the gravity acceleration, T the absolute tem- 
perature, J’ the dry adiabatic lapse rate, v the horizontal 
wind, and z altitude, measured upwards) is widely used 
in the theory of turbulence. 

2 Bannon’s definition of the horizontal shear differs 
from the ordinary one, but this fact does not change 
the qualitative reasoning. 

3 The fact that this front acts as a cold front in the 
troposphere does not alter the reasoning in the discussion 
of the earlier cases. In the complete report we shall study 
all types of fronts in connection with the proposed 
frontal model. 
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Fig. 6. Vertical cross section Aldergrove-Liverpool- 
Downham Market, 0900 GCT 25 April 1949. Heavy 
lines represent tropopauses and frontal boundaries, thin 
dashed lines are isotherms (° C), and thin solid lines 
isovels for observed wind speed (projected on a direc- 
tion of 190°, m sec). Intermediate values of tempera- 
ture are given as dotted lines. Hatched areas indicate the 
regions where clear air turbulence was observed by 
special search aircraft. 


The maximum horizontal shear through the 
front is about 60 m sec! (too km)-!, about 
the same as in the case given in fig. 3. 

The regions within which clear air turbulence 
was reported by an aircraft belonging to the 
British project, are indicated by hatched areas 
on the cross section in fig. 6. The right posi- 
tion vis-a-vis the front has been obtained 
through an estimate of the velocity of propa- 
gation of the front, obtained from consecu- 
tive sections and upper air charts. The areas 
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covered by the hatching are so km wide and 
about 300 m deep, and may be considered to 
give the true position of the turbulent area 
relative to the front. 

It is obvious that the lowest area coincides 
almost completely with the region of maxi- 
mum horizontal wind shear. 

The depth of the turbulent layers is a very 
interesting question. If turbulence were en- 
countered in connection with the frontal model 
given in this paper, it is reasonable to assume 
that the deepest turbulent layers would be 
observed at the height where the front changes 
character (in fig. 6 at about 400 mb). A careful 
investigation whether this really is the case 
cannot, however, be done by the present 
author, because of the fact that he does not 
have the material at his disposal. 

Considering the importance of the horizontal 
contra the vertical shear of the wind, the fol- 
lowing statement can be made. In Polar fronts 
as those in figs. 3 and 6, the isovels in the cross 
sections are almost parallel to the frontal bound- 
aries. This is true in the great majority of 
cases as synoptic investigations show (see e.g. 
PALMEN and NEWTON, 1950, fig. 3). 

Assuming the thermal wind equation to be 
valid in the following form (see Hewson and 
LONGLEY, 1944, p. 99—101) 


OT 
VE 5% (2) 


(where the symbols have the same meaning as 
in formula 1), we get a simple figure for the 
vertical wind shear. If we put g = 9.81 m sec~?, 


{ = 10% secrl 2>2, = 1 km, er 

Ox 
(100 km)-!, and T = 250° K, v — 4 m sec-1 
km-1.Thus, for every °C per 100 km of the 
horizontal temperature gradient the geostro- 
phic wind increases with height with 4 m sec-1 
km”. 

Assuming a front of 1 km vertical depth and 
with an inclination of 1/100, we get with the 
above-mentioned assumptions a horizontal 
(thermal) wind shear of about 4 m sec“! (100 
km)“. Thus, the often observed temperature 
difference 10°C through the front gives a 
horizontal (thermal) shear of about 40 m sec-1 
(100 km)-1, through the sloping front. This is 
in agreement with synoptically observed values 
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of this quantity, if we study the upper part of 
the troposphere. In the above discussion we 
have disregarded the conditions in the lower 
part of the troposphere. The real horizontal 
wind shear is slightly modified by the wind 
distribution in the lower layers, but the given 
order of magnitude of the wind shear is correct. 
Here it must be emphasized that with the net- 
work of wind-observing stations of to-day, the 
distance between the stations being at best 
about three times the horizontal width of the 
front, it is dangerous to assume linear distribu- 
tion of wind velocity. : 

If the above-mentioned assumptions are 
fulfilled, it is clear that strong vertical wind 
shear is connected with strong shear in the 
horizontal in the case of a normally sloping 
front, if the isovels are (nearly) parallel to the 
frontal boundaries. 

The vertical wind shear vanishes in the region 
where the temperature difference between the 
two air masses disappears, in fig. 6 at about 360 
mb, and where consequently the front becomes 
vertical. 

Before giving a short discussion of the sta- 
tistical treatment of the reports of clear air 
turbulence we want to make two assumptions. 

Firstly we assume that a great part of the 
cases with clear air gusts is connected with more 
or less distinct fronts, interpreted as in this 
paper, and the reported cases published up till 
now give the reader this impression. Having 
no opportunity to use the detailed reports from 
the search flights, we must consider this con- 
nection between fronts and clear air turbulence 
as a working hypothesis. 

In connection with this assumption a remark 
may be made concerning fig. 10 in Hislop’s 
paper. This figure gives an exemple of a jet 
stream over western Europe, and the region 
where clear air turbulence was observed is 
indicated. From the figure the reader gets the 
impression that the turbulent area is connected 
with anticyclonic wind shear. The fact is, how- 
ever, that the chart is a copy of the 300 mb 
chart, given in the London Daily Weather Re- 
ports, Upper Air Section, for 0300 GCT. The 
turbulence was observed at 1300 GCT accord- 
ing to table I in Hislop’s paper. Furthermore 
the mean height of the 300 mb surface over 
the British Isles was 8760 m, whereas the 
turbulence was encountered between 5 and 
8.6 km. As a matter of fact cyclonic wind shear, 
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and not anticyclonic, existed in the critical region at 
about 1300 GCT. 

Secondly we assume that the frontal model 
given in this paper is not too far from truth, 
and this opinion seams to be supported by the 
examples put forth in this paper as well as in 
the complete report. 

If these two conditions are fulfilled, a samp- 
ling of gust cases from all levels between 5 km 
and 12 km cannot be examined statistically 
without a special technique. Any crude method 
will certainly give disappointing results. This 
is obvious from what is said above and from 
the temperature and wind distributions given 
in figs. 3, 5 and 6. The correlation between 
horizontal wind shear, vertical wind shear, 
horizontal and vertical temperature gradient 
on one hand and clear air turbulence on the 
other cannot be studied when all reports are 
taken by the bulk. 

If clear air turbulence were reported above 
700 mb in connection with the front in fig. 3, 
one should distinguish between four different 
cases. (1) Between 700 and 500 mb the (cy- 
clonic) horizontal and vertical wind shear are 
rather considerable. (2) In the layer at about 


the soo mb level the inclination of the front 
is small and consequently the horizontal wind 
shear is also small but the vertical shear great. 
(3) A report from the region where the front 
is vertical (as in fig. 6, lower hatched area) will 
correspond to zero vertical wind shear and 
maximal horizontal (cyclonic) wind shear 
(4) Higher up the vertical wind shear is in- 
creasing but with the opposite sign compared 
with tropospheric values. 

Thus, it seems likely that any theory which 
tries to explain the occurrence of clear air 
turbulence must take all these different cases in 
account. The purpose of the latter part of this 
paper is to give a hint to a possible way of 
attacking the problem and to give some values 
of the involved parameters. 

Finally it should be pointed out that we 
might here have a possible way of attacking 
the prognostic problem of the clear air tur- 
bulence. 
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Abstract 


In earlier calculations it seems not to have been considered that the penetrating y radiation 
from the ground can, independent of the density and composition of the ground, be ex- 
pressed by the simple formula I= C x S, where S is the proportion of radioactive element 
in g/g and C is a constant. Formulas for the ionization in the case when the ground is 
covered with an absorber of greater density than air, or when the radiation travels through 
air for long distances, are derived. The ionization from various minerals is calculated and 
is found to agree in magnitude with measured values. Curves, showing the decrease in 
the ionization as a function of the depth of snow cover, are given for three different values 
of snow density. The decrease in the dose rate with the altitude above the earth’s surface 


is shown graphically. 


Introduction 


Natural ionization over the earth’s surface 
is caused mainly by radioactive radiation com- 
ing from the ground and by cosmic radiation. 
The latter radiations result in a nearly constant 
ionization, if variations during short intervals 
and changes with the altitude are not taken 
into account. In contrast, radiation from the 
ground varies widely with local conditions. 

In this communication formulas for the 
ionization due to radioactive minerals have 
been derived, and the absorption caused by 
snow and air has been taken into account. All 
the previous calculations, which were made 
30—40 years ago, were characterized by the 
fact that one (or in some cases 2—3) (KOHL- 
RAUSCH 1917) absorption coefficient was used 
to describe the total radiation from a decay 
series. A few years ago these older formulas 
were used by V. F. Hess (1947) and re- 
cently, formulas for the natural background 
radiation based on the same simplifications 
have appeared in the literature (PErRsSON and 
FRANKLIN 1951). 


Radiations 


Gamma radiation: Only the y radiation from 
the U-Ra- and Th-series and from K“° gives 
a measurable contribution to the natural ioni- 
zation, and this contribution is due to a 
comparatively small number of the y ener- 
gies. Table I gives the energy levels whose 
rhm/curie! values are greater than 0.001 and 
the corresponding absorption coefficients for 
snow and air. 

Column 3 shows the number of quanta 
emitted per primary disintegration for each 
of the energy values. This number is found 
from the decay scheme of the nuklides. For 
MsTh 2 this scheme is not known in detail. 
But about 90 % of the radiation, measured 
photographically (BLACK 1924) is due to the 
three lowest energies. If we take an average 
of these as the radiation energy of MsTh 2, 
the error in the final value will be negligible 
since the radiation dose from MsTh 2 amounts 


' rhm is the radiation in r units per hour at 1 m distance, 
r being the international dose-unit measured by ionization 
in air. 


to about 5 % of the total dose from the tho- 
rium series, which means approximately 1% 
of the total natural radiation on the earth’s 
surface from the most common minerals (see 
below). 

Several schemes for the disintegration of 
K*° have recently been proposed (SPIERS 1950, 
MATTAUCH and FLAMMERSFELD 1949, GRAF 
1951). Since none of these seem to have any 
marked superiority over the others, the calcul- 
ations have been made for three different alter- 
natives. 

Actinium gives no measurable radiation 
because of its very small occurrence and the 
low gamma energies. C™, and probably Rb:7 
are beta emitters. 

Beta radiation: The B radiation is absorbed 
by the air already at a few meters above the 
earth’s surface. The ionization due to f radia- 
tion is usually small because only those particles 
emitted from a thin layer near the surface 
actually escape into the atmosphere. Even 
with very thin-walled measuring chambers the 
ionization due to ß rays from the earth’s 
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Table I 
mn i ee 
n & 
hv uanta hyn i = Mr : 
Element Mev a pri ee snow snow snow ol 
disintegr. Q = 0.10 — 02 0 = 0.40 
Ra 0.184 0.012 0.0005 1.240 10 | 3113 10021 41080. 1702|, 16.12. 1025 
RaB 0.241 0.1IS 0.0066 1.150 2.873 4.603 14.88 
0.294 0.258 0.0205 1.078 2.692 4.314 13.94 
0.350 0.450 0.0457 1.008 2.519 4.037 13.04 
RaC 0.607 0.658 0.1454 0.801 2.000 3.205 10.36 
0.765 0.065 0.0199 O2 1.805 2.892 9.346 
0.933 0.067 0.0270 0.659 1.647 2.639 8.528 
1.120 0.206 0.1061 0.599 1.496 2.398 7.749 
1238 0.063 0.0373 0.563 1.406 22265 7.282 
1370 0.064 0.0433 04527 1.316 2.109 6.815 
1.761 0.258 0.2362 0.439 1.098 1.759 5.682 
2.198 0.074 0.0887 0.376 0.940 1.506 4.868 
MsTh 2 0.090 I.00 0.0113 1.514 3.783 6.061 19.59 
ThB O.IISI 0.03 0.0005 1.451 3.625 5.808 18.76 
0.2381 0.91 0.0516 1.156 2.888 4.627 14.95 
0.3001 0.03 0.0024 1.069 2.670 4.278 13.82 
The” 0.2518 0.02 0.0005 1.129 2.820 4.519 14.60 
0.2767 0.10 0.0025 1.096 DTA 4.386 TAT, 
0.5110 0.17 0.0103 0.861 2.051 3.446 IT.I4 
0.5832 0.60 0.0441 0.816 2.038 3.266 10.55 
2.624 0.92 0.4782 0.328 0.820 7313 4.244 
3.240 0.08 0.0540 0.271 0.677 1.085 3.505 
Kee 1.54 0.04 0.0310 0.488 1.218 1.952 6.308 


surface will be very small and will not be 
considered here. 

“Bremsstrahlung” : However, ß radiation gives 
a secondary effect, the “Bremsstrahlung”, 
which must be taken in account. The approxi- 
mate formula 


(dE/dx) “Bremsstrahlung”’ 


(dE/dx) ionization 


NET 


820 


shows that for K4° in the most common 
igneous rocks, for which the average of the 
atomic number is about 12, approximately 
10 % of the energy of the ß radiation is trans- 
formed into the “Bremsstrahlung’’. Its total 
absorption coefficient is somewhat higher than 
that of the y radiation. The “Bremsstrahlung” 
dose from K*° on the earth’s surface may be 
about 10 % of the y dose from the same isotope. 

A rough estimate of the contribution of the 
“Bremsstrahlung” to the ionization due to 
radium and thorium in the ground shows this 
to be 1 % or less of the total ionization from 
these elements. 
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The absorption coefficient 


The total absorption coefficient y is given by 

u =t+o-+2, where 7 is the absorption 

coefficient due to the photo effect, o to Comp- 

ton scattering and absorption and x to pair 

production. The energies recorded in table I 

vary between 90 and 3,240 keV. Between 200 

and 2,000 keV, energy absorption due to 

photo effect and pair production is negli- 

.gible compared with that due to the Compton 

effect for the absorbers considered below. 


The photo effect is of secondary importance 
compared with the Compton effect even at 
energy levels as low as 90 keV (MsTh 2) in the 
absorbers considered. Since the radiation from 
MsTh 2 contributes only about 1% of the 
total absorbed energy it is not necessary to 
take the photo effect into account. 


Pair production, which occurs for energies 
greater than 1.02 MeV, is measurable only for 
the two highest y energies of 2.624 and 
3.240 MeV respectively from ThC”. For these 
energies x is 4 % and 7.4 % respectively of o. 
Taking these corrections into consideration a 
value for the dose rate is obtained that is 
about 99.7 % of that without correction. 
Consequently, the pair production can also 


be neglected. 


The Compton absorption coefficient can be 
calculated by the use of KLEIN-NISHINA’s (1929) 
formulas. The total Compton coefficient per 
electron, ,0, and the Compton absorption 
coefficient per electron, .o,, are given as 
functions of the energy in MeV in fig. 1. 


The third curve of fig. r shows the depend- 


o 
ance of (22) = f; on the energy. To find 
el Ji i 7 


the absorption coefficient we have to multiply 
the Klein-Nishina coefficient with the number 
of electrons per cc. We arrive at 


Zi 
Oj == L 0,0 Zp: = 


vr 
W, and o, = Lee pi Tr 


W; 


where L is Avogadro's number (6.023 + 1024) 
o the density 


Z; the atomic number of the element 


with index i, contained in the 
absorber 

W; the atomic weight of the same 
element 
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Fig. 1. The total Compton coefficient per electron, 

«0, the Compton absorption coefficient, eda, and their 

quotient eOa/eO are given as functions of the energy 
in MeV. 


pi the part by weight of the absorber 
that consists of the element with 
the atomic number Z; 


Scattered radiation 


The lack of information on the scattered 
radiation makes calculations of this kind 
approximate. This problem has been dealt 
with in a theoretical way by several authors 
(HIRSCHFELDER et al. 1948, Lewis 1950, MER- 
EDITH 1951), with different limiting assump- 
tions. Using values published by HirscHFELDER, 
MAGEE and HULL (1948) for the part of the 
total intensity which is due to unscattered 
radiation, a rough estimate of the contribution 
of the scattered radiation has been made in the 
following way. Let us imagine a series of 
planes in the ground parallel to the earth’s 
surface separated by a distance of one cm. 
The ionization due to the radiation from a 
layer bounded by two of these planes can be 
expressed in percent of the total primary radia- 
tion. From the values given by Hirschfelder et al. 
we find what part this is of the total radiation 
from the layer in question. Dividing by this 
number we get the total radiation from the 
layer. By adding the radiation from the 
different layers an approximate integration is 
made. The total radiation determined in this 
way will be expressed in per cent of the primary 
radiation. These calculations have been made 
for 1 MeV and for 3 MeV. In the first instance 
the total radiation was found to be 118 % of 
the primary one and in the latter 116 %. 
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For this estimate the following simplifica- 
tions have been made. Hirschfelder’s formulas 
are valid only for a homogeneous bundle. As 
will be seen below, the natural radiation from 
the ground in special cases behaves like a 
homogeneous bundle, for if the absorption 
caused by the air is neglected, the dose is not 
dependent on the distance from the earth’s 
surface. This simplification makes the values 
too small. However, no account is taken 
of the fact that the scattered radiation is 
softer than the unscattered one, and this 
will make the correction too high. Since, 
as was said above, the effect of the scattered 
radiation is almost the same, expressed in 
per cent, for radiation of the energies one and 
three MeV, and because the energies of greatest 
importance in the practical case are of a magni- 
tude between these limits, we see that we can 
take the correction factor to be 1.2 and that 
we can be almost sure that this figure should 
not be I.I or 1.3. 

A correction of 20 % for the scattered radia- 
tion is therefore added to all values obtained. 
The calculated radiation for potassium is 
corrected by an additional 5 or 10 % to 
account for the “Bremsstrahlung”’. 


Theoretical calculations 


Fig: 2. 


The ionization contribution in point 
P (fig. 2 above) from the ring-shaped 
volume element dV and of energy E; is 


I Ejn; Ojai 


dl; = cdV - 
rer. 96.77. 10% 
. 2e A 1010 “ e-Gig di + Fj) dz) 
roentgen/sec. 
I; = dose rate due to quanta of energy 
E; 
r = distance in cm 
E; = quantum energy in MeV 
n; = number of quanta of energy E; 


per primary disintegration 

ia = absorption part of the Compton 
coefficient for the air and energy E; 

oi = total Compton coefficient for the 
air and for energy E; 

dig — total Compton coefficient for the 
ground and energy E; 

c = number of curies of radioactive 
element per cc 


Here 3.7 : 101%. c is the number of disinte- 
grations/cc sec. Let us assume that the ground 
contains S g radioactive element per g. The 
amount of the nuklide per cc is Seg, where 0, 
is the density of the ground. The number of 
atoms of the nuklide per cc is Se,L/W where 
W is the atomic weight and L Avogadro’s 
number. The number of disintegrations per 
second will be 
a. SaL _ 09.693 . 80,1 

W Tr W 


= 37 RO TEE 


Inserting this expression for c, adding over 
all energies and multiplying by 2.08 - 10°, 
we find the ionization in ion pairs/cc sec.: 

5% I 


Tiers Hate 
- Z'E;n;oia * ehrt od). AV 
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ion pairs/cc sec. 


Here is 


dV =2nbdbdz, d = - 


- r and 
== 


The radiation over an infinite plane is expressed 
by the following, if we take into account the 
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radiation coming from layers down to a depth 
of A cm under the earth’s surface. 


d en a—z dé 
seman fee [5 mes 15e 
If we take r as the independent variable and 


take 


AO — ZOig 


- r= y, the expression becomes 
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where Ei denotes the exponential-integral. 
For this is valid 


Sera 1414 
2 d = =e mb 
t b+ 


eb 
—= — Eil— er + b)] 


== —22,030..10°°-» 


er mt 


Thus 
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where ¢ (x) = fz is 


x 


By partial integration we find 


y . 
J Eil— m(b + 1))d = — e—™ [x - d(x)” + 
Lei DR 
ee ‘fx À ie == (y + 4) + 
- Eil— m(y + b)| —b - Ei(— mb) — 
pe etalon ie ey 


m 


As y approaches infinity this becomes 


co e mb 


f Eil-m(b + 1)]d= —b- Ei(— mb) — — 


o m 


or in the practical case 
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If the absorption caused by the air is neg- 
lected, which is permissible if the distance 
from the earth’s surface does not exceed a 
few meters, we get 


S Og Ojal 


102%. S SIE ERS 
Tyears W i Gig 


Bor On 


ion pairs/cc sec. 


We see from this that the dose rate is inde- 
. > 
pendent of the distance from the earth’s sur- 
face, taking into consideration the limitations 
imposed by the absorption due to air. 
If we introduce the expressions for oj, and 
Sia from page 56 this becomes 


So; 
f=2.07"= 167% = 
Es W 
FR 
Z pil a 
J i dE; n;f; 
Ep, Zug i 
y = Wy, 


0 
where f; = (<2 


a . . 
| of the energy in question 
e 1 


(see fig. 1) and o, is the density of the air. For 
the “elements contained amie ‘and most com- 


; i 7; 
mon minerals, Z'p - W has a value near to 


0.500. For the minerals for which calcula- 
tions have been made the value is between the 
limits 0.497 and 0.500. Therefore the formula 
can be written 


Minis 
vets W 


ion pairs/cc sec. 
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Fig. 3. The upper curve refers to the upper scale on 


the abscissa and the lower curve to the lower one. 


Z'E;n;f; is a constant. Thus the radiation from 


t 
a decay series can, near the earth’s surface, be 
expressed by I = const. - S 

It is apparent that the only variable is the 
proportion of the radioactive element ex- 
pressed in g/g. We need know neither the 
density of the ground nor its composition 
(and so the absorption coefticients) to be able 
to calculate the ionization on the surface. 

If, however, the absorption caused by the air 
cannot be neglected or if the earth’s surface is 
covered with an absorber of greater density, 
for example snow, the formula must be written 


S 


I = 2.63 - 1017 + ———— 
: Fa 


» SE ni figi 


where 


g = ao; Ei (—a0;) + ei 


or if we write ao; = x; g;= x - Ei(—x) +e* 
where o; is the total Compton coefficient of 
the absorber. g (x) is seen in fig. 3. The upper 
curve refers to the upper scale on the abscissa 
and the lower curve to the lower one. 

If we takc’into account the absorption caused 
both by the snow and the air g; has the form 


gi = [bo; + (a — bo] - Ei[— bo; — 


= (a =— b) oil] Sr e— 49;—(a—B) 951 


(where b is the depth of the snow cover and a 
is the distance from the earth’s surface to 
the measuring point) or by writing bo; + 
+ (a— boy = x 

gi NT Ei (— x) + Cu 


This function is also valid if the air 
absorption it neglected, and we see that 


& (cs) gl) = g (1 + Xo). 


Practical calculations 


Because the statements as to the half-life and 
decay scheme of K4° differ so widely, dose rate 
calculations for the potassium radiation have 
been made for three alternatives. 

The expressions for the dose rates due to the 
different decay series are: 


1 Os3O5 ° Moto Sn, 


Im= 0.534 10° + Sty 
a. Ix = 44.1 . Qype 
b. Ik = 27.3 SRE 
G Ike = 30.77 I 


The alternatives a, b and c refer to the follow- 
ing constants: 
d 12.0.5510 hyearsyy Wera 
bovis =) 1. treo years Gu = 5? 
1.33%. 107 VCATS Vo ei 
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Fig. 4. The decrease in the ionization with the depth 
of the snow cover is given for three different densities 
of the snow. 


1 Sprers (1950). ? MATTAUCH and FLAMMERSFELD 


(1949). © GRAF (1951). 
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The potassium radiation is corrected by 10 % 
for the “Bremsstrahlung” in the alternatives a 
and b and by 5 % in the alternative c. It is 
assumed that 0.0119 % of all potassium is K*° 
(NIER 1950). 

The proportions of radium, thorium and 
potassium (observe that it is the usual potassium 
isotope) in some common minerals are collected 
in table II. 

On the basis of these proportions the values 
of ionization over different minerals in table III 
have been calculated. 

Fig. 4 shows the decrease in the ionization as 
a function of the depth of the snow cover, for 
three different values of the density of the snow. 

It is to be observed that the curves show only 
the decrease in the radiation from the ground. 
The almost constant contribution of the cosmic 
radiation makes the observed variations less. 
The earth’s surface is assumed to be composed 
of igneous rocks having an average composi- 
tion. The form of the curves in the case when 
the ground has another composition are very 
similar to that of fig. 4. 

The decrease in the radiation with the alti- 
tude above the earth’s surface is given in fig. 5. 
(The earth’s surface is assumed to be composed 
of igneous rocks.) The lower curve shows the 
decrease in the primary radiation. The upper 
one gives an upper limit for the radiation, 
taking into account only the absorption part 
of the Compton coefficient of air in the ex- 
ponential factor. This means that it has been 
assumed that all scattered radiation is directed 
forward and is of the same wavelength as the 
unscattered radiation. 

All values required for the calculations are 
given in table I. Sprers’ values for the con- 
stants of K4° are used in the calculation of the 
values on which the curves of figs 4 and § are 


based. 


Discussion 


The errors depending on the uncertainty of 
the constants used are only 2 to 3 % for Ra 
and Th, whereas they are 40—so % for K40. 
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Fig. 5. The lower curve shows the decrease in the 


primary radiation with the altitude above the earth’s 
surface. The upper curve gives an upper limit for the 
radiation. 


This is evident from the expressions for the 
three alternatives for the potassium radiation. 

If we assume that one of the three statements 
as to the constants of K*° is correct, we can say 
that the error of that value of the total dose 
rate which is obtained by use of the corre- 
sponding expression for the radiation from K®®, 
after correction for scattering and “Brems- 
strahlung”, is 15 % or less. This is the case if, 
as we can assume, the error of the correction 
factor for the scattered radiation is less than 
10 %, the error of the primary dose is 2—3 %, 
and if that of the correction factor for the 
“Bremsstrahlung” from potassium results in 
an error of 1—2 %, in the final value. To this 
are to be added the errors caused by the un- 
certainty of the proportions of radioactive 
elements given in table II. 

In a recent paper some measurements have 
been published (SIEVERT and Hurtovist 1952) 
showing a variation in the radiation over 
different rocks from 5.5 to 10.5 ion pairs/cc 
sec., and these values are situated between the 
limits which are to be expected from table II. 
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Table II 


Mineral 
pens) ero. per g) - 108% 


(g Ra226 | (g Th?23?2 | (g K39 


Igneous rocks: 
TASTES N CU uen ee 143 12 2.6 


Acidic rocks: > 65 % SiO, 


GRAMME CT A DE ANR EN EEE LE RE. 3.1 20 3.4 

QUney ranted A ee Rois es 1.0 9 3.9 

Younger granites in Finland (max. values)? .... 6.5 59 Si 

Grenodiorit? car, NO ens «ke heen Mic bars 227 18 2.5 
Medium rocks 65—55 % SiO, 

[BEGGAR LOTS Gres EN RM CARN RO 1.4 6 1.7 
Bastenrocksı< ss) LSIO, 

CADÈTO PEN er cyt. Le SR MER. 2 0.83 5.1 0.7 
Ultrabasic rocks 

Bene HEC RIO es OLR Sa EEE 0.52 3.3 0.8 
Sedementary rocks : 

IS a NRA ORNE, Fe AR sie. 1S I = * 

SAROS ONE RE A ee tee ee EH TA > | Sue ID 

ERNEStCONSEE IRs Big ice ehr 0.94 RE - OB 


1 RANKAMA and SAHAMA 1950. ? HOLMES 1930. 3 Hess 1947. * HOLMES 1915. 5 MACHATSCHKI 1948. 


Table III 
Calculated values of the dose rate over various minerals 


The dose rate is expressed in ion pairs/cc sec. For each mineral the values corrected for scattered radiation 
and ‘‘Bremsstrahlung”” are given by the numerals in bold face. The cosmic radiation is not included. 


Ionization due to content of Total ionization 
Mineral 
Ra Th | Ka K? | Ke 1 | 2 3 
Igneous rocks: 
ERR STEELS SS pet TR CO CRO OUI CREO NEE 07, 0.6 Te 027 13 2.5 27 2.6 
: 3.1 .6 3.2 
Acidic rocks: > 65 % SiO, 
(GRAN EN eco 108 Tall re 0.9 1.7 4.4 3.9 4.7 
| 5.4 4.8 5-7 
ONE Water, AE PR EE ere te 0.6 0.5 7 1.0 2.0 2.8 2.2 GR 
3.5 2.7 3.8 
Younger granites in Finland (max. values),.| 3.7 32 262) 1.4 2.6 9.2 8.4 9.5 
1124 770.22 Urs 
Gran oO done ee are cc IS 1.0 DT 0.7 3 3.7 273 3.9 
PR 1 4.016748 


Medium rocks 65—55 % SiO, 
TD YG YB ey Ce mean ices CEG mous SOS sa ach 0.8 0,32 07 0.46 0.9 


Basic rocks <,.55% SiO; 


CAD OM Mer A ee labs ee 0.5 0.27 0.31 0.19 0.35 ig! 1.0 OR 
1.2 1.2 13 
Ultrabasic rocks 
PETIAOGIt ee ee tod eh wie nie à 20 0.30 0.18 0.35 0.22 0.41 0.8 0.7 0.9 
1.0 0.8 Let 
Sedementary rocks: 
SAandStONE Le ce cle mc es See ses atte ara aan 0.8 0902 12.049. |20,302 11.0.0 1.6 1.4 1.7 
1.9 ne? Pe 
TRES CON CR een ie Cheney ecu ce 0.5 0.05 0.13 0.08 0.15 0.7 0.6 0.8 
0.8 0.7 1.0 
CHERE TRS A ro OR, 0.9 0.7 Tek 0.7 13 207 270) 2.8 
3-3 2.7 3.4 
EEE Be EE EEE ee ee ee 


1 Spiers 1950. 2 MATTAUCH and FLAMMERSFELD 1949. © GRAF 1951. 
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Experiments on Non-Laminar Flow of Mercury in Presence 
of a Magnetic Field 
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Abstract 


Two examples are given of non-laminar motion in an electrically conducting liquid in 


the presence of an external magnetic field. 


Introduction 


In many cosmical applications, e. g. in the 
interior of the sun, in the interstellar space, 
in the interior of the stars, in the ionosphere 
and when cosmic radiation is concerned, the 
presence of a magnetic field gives rise to forces, 
which are comparable with other forces. The 
purpose of the following discussions is to give 
some special examples of these effects, acting 
upon a liquid which is in a state of non- 
laminar motion. 

The behaviour of an incompressible liquid 
conductor is strongly influenced by the pres- 
ence of an external magnetic field. This is 
due to the fact that every motion perpendicular 
to the magnetic field will induce electric 
currents which interact with the magnetic 
field, resulting in forces on the elements of 
the liquid. In this way it can be shown that 
one layer of the liquid is coupled to the 
neighbouring layers, not only by viscous forces 
but also by electro-magnetic forces. 

The laminar state of motion in such a 
liquid has been treated theoretically and ex- 
perimentally by HARTMANN (1937), ALFVÉN 
(1942), WALÉN (1944), LUNDQUIST (1949) and 
many other authors. On the other hand, 
until now, little attention has been paid to the 
treatment of a non-laminar state. HARTMANN 


(1937) gives some aspects on this subject and 
the spectrum of magneto-hydrodynamic tur- 
bulence is discussed by BATCHELOR (1950) and 
CHANDRASEKHAR (1951). 

While studying the flow of mercury in 
channels placed in a strong magnetic field 
HARTMANN (1937) found that the transition 
between laminar and non-laminar motion as 
well as the non-laminar state itself were in- 
fluenced by the field. The measurements with 
rotating cylinders, described in Sec. 1, are in 
agreement with Hartmann’s results as regards 
the non-laminar state and the explanation 
seems to be the same as pointed out by Hart- 
mann, viz. a suppression of the whirls by the 
magnetic field. A more detailed discussion is 
given in a previous paper (LEHNERT, 1951). 

Sec. 2 gives an example of a vortex-configu- 
ration which in some respects may be similar 
to KARMAN’S vortex-street (KARMAN, IQIT). 


1. Measurement of the Apparent Viscosity 
with Rotating Cylinders 


The measurement was carried out with the 
well-known method of balancing the torque 
between two concentric, non-conductive cyl- 
inders, one of which is at rest while the other 
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Fig. 1. Co-axial cylinders rotating with relative angular 
velocity wo. Magnetic field, Bo, in axial direction, and 
particle-velocity, 7, in a plane perpendicular to the axis. 


is rotating with constant angular velocity wo. 
The space between the cylinders was filled 
with mercury and the external magnetic 


field, Bo, parallel to the axis (fig. 1). The 
torque acting on the cylinder at rest was indi- 
cated by an elastic wire and a mirror, reflecting 
a light-ray on a scale. fig. 2 shows the result 
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Fig. 2. Torque in newtonmeters + 10° (= dyn.cm) as a 
function of angular velocity. Outer cylinder rotating. 
1; = 37.0 mm; ro = 50.5 mm; I, = 9 mm; l,= 78 mm. 
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with the outer cylinder rotating and fig. 3 
with the inner cylinder rotating. : 
The curves for the purely mechanical case 
differ distinctly from each other. According 
to TAYLOR (1922) the instability arises earlier 
when the inner cylinder is rotating than in the 
opposite case. It may be convenient to distin- 
guish between two domains of the curves: 


(i) Small velocities. The presence of a 
magnetic field gives rise to electric polari- 
zation. The velocity-distribution is prac- 
tically undisturbed in the upper regions of 


the cylinders, where the velocity v varies 
insignificantly with the coordinate in the 
direction of the field. In the region 
between the lower edge of the inner cyl- 
inder and the bottom, however, the distri- 
bution of the mechanical case is strongly 
modified by the field resulting in a great 
increase in torque, as long as the motion 
is laminar. But if this ‘“edge-effect” alters 
the velocity-distribution a non-laminar, 
thin boundary-layer may arise at the walls 
in presence of the field, even if the mechan- 
ical case is laminar. This may explain 
the fact that the increase in torque in 
reality becomes rather moderate. 


(ii) Large velocities. The edge-effect will 
increase the torque both in a laminar and 
in a non-laminar case. But if the motion is 
non-laminar the field may have a second 
effect in respect to the torque, acting in 
the same or in the opposite direction as 
the edge-effect. If it acts in the same direc- 
tion the field will always increase the 


Torgue 


2 


RR 0.32 Wem? I20090055 


600 | 07) 


500 | 30) 


300 
| 2 
300 
200) 
Anguior velocity 
(‘sec’) 
sie iG 5 7] 75 20 + 


Fig. 3. Inner cylinder rotating. 
rj = 37-0 mm; ro = $0.5 mm; |, = 11 mm; /, = 95 mm. 
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torque. If it acts in the opposite direction, 
however, the net effect will be an increase 
or decrease due to the dominating one of 
the effects above. fig. 3 shows that the 
net effect is a decrease, which makes the 
assumption very plausible that a non-laminar 
motion of a conducting liquid is suppressed by 
an external magnetic field, at least in certain 
cases. 


2. A Magneto-Hydrodynamic Vortex- 
Configuration 


If a motion is started in a liquid with infinite 
conductivity and no mechanical viscosity, the 
presence of an external magnetic field will 
force the motion to be propagated along the 
field-lines leaving the regions outside undis- 
turbed. In a liquid with finite, but large, con- 
ductivity and small mechanical viscosity the 
motion will spread outside the fluxtubes which 
contain the region where the motion is started, 
but if the fluxtubes of the liquid considered 
are not too long the behaviour is roughly the 
same as in the former case. 

A special form of motion is studied with 
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Fig: 4. Arrangement for producing the configurations of 
Figs. 5 and 6. The strap is driven by a motor. Field- 


strength Bo & I Vs/m? (= 104 gauss). 


$— 200579 


——— 40mm ee 


nm 60mm 


Fig. $ a. Outline of the surface as observed when looking 
from the top of the vessel in fig. 4. The dark region 
follows the motion of the copper-disc. The small ar- 
rows indicate the direction of rotation of the whirls. 


Fig. 5 b. Cross-section of the surface. The dark region 
corresponds to the layers following the motion of the 
copper-disc. The surface is horizontal in the regions 
at rest and the whirls are indicated by the cavities at the 
boundaries. 


the apparatus of fig. 4. At the bottom of the 
vessel a rotating copper-disc, driven by a 
motor, is situated. The surface of mercury is 
about 15 mm above the copper-disc. Imme- 
diately above the latter a thin, immovable 
disc of non-conducting material is situated, 
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Fig. 6a. Picture of the surface when the disc is being 
accelerated. 

b. Stationary motion, @o about 10 sec.~. 

c. Disc suddenly stopped. In this case five large whirls 
moved towards the centre. 
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which screens off the inner part of the copper- 
disc. The magnetic field is normal to the discs. 

If the magnetic field is absent no mo- 
tion is observed at the surface when the 
copper-disc is set into rotation. With a mag- 
netic field-strength of 1 Vs/m? (= 104 gauss) a 
motion is observed on the surface as shown 


by fig: 4: 


(i) The region outside the radius of the 
copper-disc as well as the region which is 
screened by the immovable disc are at 
rest, whereas all layers situated between the 
radii of-the two discs form a cylindrical 
sheet, moving with the same angular 
velocity as the copper-disc. The dark 
ring in fig. 6 b corresponds to this region. 


(ii) If the linear dimensions are suitably 
chosen two vortex-rows are superposed 
on the motion of the cylindrical sheet and 
the centres of the whirls are situated above 
the edges of the discs. The rows move 
slowly along the sheet and the rotation is 
directed in such a way that they seem to 
roll like a ball-bearing between the moving 
and immovable regions. 


Fig. 6a shows the surface when the disc is 
being accelerated and the whirls appear di- 
stinctly. In fig. 6b the motion is stationary 
with a rather high angular velocity (about 
10 sec.”!) and fig. 6c shows the motion of 
some large whirls towards the centre as the 
disc is suddenly stopped. Stationary rotation 
at small angular velocities gives a picture 
similar to fig. 6a. Ata first glance the con- 
figuration seems to be similar to Kärmän’s 
vortex-street (KARMAN 1911) in pure hydro- 
dynamics. Karman has shown that the only 
stable configuration of vortex-rows in a purely 
hydrodynamic case consists of two rows, the 
centres of the whirls in one row being situated 
in front of the spacing between the whirls 
in the other. There are, however, at least two 
differences between Kärmän’s case and that of 
fig. 6: 


(i) The positions of the whirls in one row 
seems to be independent of the others as 
shown by fig. 6a. 


(i) A stable system consisting of only one 
vortex-row exists in the magneto-hydro- 
dynamic case; if the immovable disc is 
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taken away the outer row still exists, 
which is distinctly shown by using a 
copper-disc of about the double radius of 


that in the experiment described. 


Both cases above are clear examples of 
the importance of a magnetic field for the 
non-laminar motion in media with as good 
electrical conductivity as mercury. In cosmical 
applications, where the linear dimensions are 


much larger, the conditions for magneto- 
hydrodynamic phenomena are much more 
favourable than can ever be realized in the 
laboratory, and the effects will have a still 
greater influence on the motion. 

The author wishes to express his thanks to 
professor H. Alfvén who has suggested this 
work, and which has been made possible by a 
grant from Statens Naturvetenskapliga Forsk- 
ningsrad. 
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Comment on a Note on the Kinetic Energy 
Balance of Zonal Wind Systems 


By J. van MIEGHEM 
University of Brussels 1 


March sth, 1952 


Since Kuo’s equation of balance of the 
kinetic energy of zonal motion in the at- 
mosphere is of basic importance for the 
study of the general circulation (Kuo 1951), 
I wish to point out that this equation may 
be simplified if account is taken of the equa- 
tion of continuity for the mean motion of 
the atmosphere. It turns out that the con- 
stancy of the specific mass o along latitude 
circles and with time appears to be an un- 
necessary hypothesis. Moreover the terms 
related to the meridional circulation are more 
easely interpreted. 

In orthogonal and axial symmetrical co- 
ordinates (x? = 2, x?, x8), the equation of 
motion along a latitude circle assumes the 
form of a balance: 
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from which the equation of the mean mo- 
tion may easely be deduced: 


where X designates the zonal mean of X 
and X the corresponding weighted mean 
(oX = 0X), X” = X—X the fluctuation 
of X with respect to the weigthed mean, 


t At present: Institute of Meteorology of the Uni- 
versity of Stockholm. 


yy the coefficients of the metric form (y;; + 0, 
yg =O when i + j;'i, j= 7 2,3) Ne 
the contravariant and v; = y;; vi the covariant 
components .of the air velocity v relative 
to the Earth, y the determinant || yj; ||; where 
Or OUR 0. 
= Spot 
dt at ou 
when i= 1 or + I. 

After multiplication of (1) by v! and 
substitution of the equation of continuity 
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or, explicitly in spherical coordinates (y,, = 
ei, Divan + PAR ss > o I 7) 
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dy r A 
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= = 5): using the classical zonal, meridi- 
OZ or 

onal and vertical components (u, v, w) of 
the velocity, 
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are the meridional and vertical components 
of the flux of the “mean zonal kinetic en- 
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is the corresponding rate of production. 
The flux (Y, Z) of “mean zonal kinetic 
energy” in the meridional plane includes a 
convective flux and a flux related to the 
flux of linear zonal momentum; this non 
convective flux is directed to the North in 
the belt of westerlies (VAN MIEGHEM, 1949). 
The production of “mean zonal kinetic 
energy” results from: a) the mean meridi- 
onal circulation (@v, ow) and the vorticity 
created by the Earth’s rotation and by the cur- 
vature of the zonal flow; b) the eddy stres- 


ses ou”v” and ou" w" and the horizontal and 
vertical shears of the mean zonal wind n. 

It should also be noted that the “mean 
zonal kinetic energy” (& (a)? = =v, n) is 
not a scalar quantity, (but disregarding the 
factor 3) one of the components of the ten- 
were (1,7 == 1, 2, 3). 

Finally Kuo’s conclusive statement “the 
eddy process in the atmosphere is a self con- 
sistent and self sufficient mechanism in main- 
taining the zonal wind systems” may be mis- 
leading and instigate unnecessary controverses. 
In order to avoid misunderstandings, it should 
be mentioned that the rates of production 
(o) of the potential (®) and kinetic (Kn) 
energy of the mean three-dimensional mo- 
tion, of the kinetic energy (K.) of the large- 
scale eddies and of the thermal energy (E) 
(sensible and latent heat) assume the form 


En + A — co(®), 
o(K,) = — A+(-Vp-v’) , 
ae =—pdivv —(-Vp ") z 


where 


For simplicity the viscosity and the small- 
scale turbulence have been neglected here. 
The above formulas show that: 


1) p divv and — VY p : v” express the 
rate of conversion of thermal energy resp. 
into kinetic energy of the mean motion and 
into kinetic energy of the large-scale eddies ; 

2) /\ expresses the rate of conversion of 
kinetic energy of the large-scale eddies into 
kinetic energy of mean motion, or, more 
precisely, a positive term of /\, such as 


A " où + u t = a t 
De eS = O, gives a rate 
Q dy À sp 3 


of conversion of K, into K,, and, a negative 


— /odii n 
term, such as ou"w" [= — —) <o, a 
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rate of conversion of K,, into K,. 


Conversion of thermal energy into kinetic 
energy of the large scale eddies occurs when the 
work performed by the pressure force along 
the path of the large-scale eddying motion 
is positive; when this work is negative 
(stable large-scale eddying motion) the con- 
version of energy takes place in the opposite 
direction. 

In a three-dimensional diverging flow, 
thermal energy is transformed into kinetic 
energy Km; the reversed transformation is 
taking place in a three-dimensional con- 
verging flow. 

From the point of view of large-scale 
energy transformations in the atmosphere, it 
should be emphazised that all possible con- 
versions between the mechanical energy 
Kn + ® of the mean motion, the thermal 
energie E, and the kinetic energy K. of 
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the large-scale eddies may occur and in fact, 
do occur. 
Stockholm, February, 1952. 
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Reply 
By H. L. KUO. 
Massachusetts Institute of Technology 


March 1952 


Prof. Van Mieghem’s introduction of the 
mean value weighted by the density distribu- 
tion is interesting. Through this definition of 
the mean, the restriction on the density put 
forward in my short note becomes unnecessary. 


As for the sentence “the eddy process in the 
atmosphere is a self consistent and self sufficient 
mechanism in maintaining the zonal wind 
systems” quoted by Prof. Van Mieghem, I 
have only intended to point out that this eddy 
process is important both for the angular 
momentum balance and for the kinetic energy 
balance of the mean zonal winds and that its. 
effects are of the right order of magnitude and 
direction, so that it may account for the major 
portions of these balances. It is hoped that it 
does not give the impression that no energy 
supply is needed by these eddies, nor that it 
excludes the effects of other physical processes. 
The note is primarily concerned with the 
kinetic energy balance of the mean zonal wind 
systems, in which the work done by the 
pressure forces does not occur, although it is a 
very important factor for the production of 
kinetic energy of almost all other kinds of 
motions. This is one advantage of my equation. 
The energy exchange between eddies of 
different scales discussed by Van Mieghem is 
surely an important process in the atmosphere, 
but not directly connected with the main- 
tenance of the mean zonal wind systems. 
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